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ABSTRACT 


Two  single-stage,  transonic  compressor  designs  were  tested 
under  various  undistorted  operating  conditions  to  characterize 
the  process  leading  up  to  aerodynamic  stall.  In  addition  to 
extensive  instrtimentation  to  measure  the  compressor  steady 
performance,  the  rig  case  was  instrumented  with  eight  high- 
response  pressure  transducers  equally  spaced  around  the  annulus 
for  stall  development  detection.  High-response  measurements  were 
low-pass  filtered  and  both  spatially  and  temporally  analyzed 
using  discrete  Fourier  techniques . 

At  all  speeds  tested  for  both  designs,  stall  inception  was 
characterized  by  growth  of  a  small  amplitude  rotating  wave.  The 
waves  did  not  grow  signifcantly  until  just  prior  to  the 
instability,  when  exponential  growth  into  fully-developed 
rotating  stall  occurred  very  rapidly,  within  6-10  rotor 
revolutions .  The  amount  of  time  the  rotating  waves  could  be 
detected  prior  to  stall  varied  considerably  with  compressor 
operating  condition  and  was  largely  dependent  on  the  local  slope 
of  the  compressor  speedline  characteristics.  Stall  warning  times 
ranged  from  less  than  one-tenth  of  a  second  to  more  than  two 
seconds  for  the  same  machine  operated  at  different  high  speeds. 
The  influence  of  compressibility  effects  are  also  discussed. 
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1 . 0  INTRODUCTION 


1 . 1  Motivation 

The  search  for  an  incipient  stall  signal  in  aircraft  gas 
turbine  engine  compressors  is  not  new.  The  interest  is  primarily 
motivated  by  a  desire  to  actively  control  compressor  aerodynamic 
instabilities  (stall) .  Because  of  the  need  to  avoid  these 
instabilities,  substantial  stall  margin  (or  stability  margin  or 
surge  margin)  is  designed  into  fan  and  compressor  components. 

The  margin  is  required  for  the  engine  to  maintain  adequate 
performance  under  the  most  severe  aircraft  engine  transients  and 
flight  conditions.  As  a  result,  component  performance  is 
sacrificed  at  more  favorable  engine  operating  conditions.  A 
reliable  indicator  of  an  impending  instability  (stall  precursor) , 
if  detected  with  sufficient  warning  time,  would  allow  for  active 
controls  to  take  action  to  avoid  or  suppress  the  stall.  In  its 
crudest  form,  this  active  stall  control  would  use  bleed  air, 
varizd>le  blades,  and/or  fuel  control  to  back  away  from  the 
insteUsility .  In  a  more  elaborate  form,  coo^ressor  \2nsteady 
performance  would  be  controlled  to  provide  a  dynamic  con^ressor 
stabilization.  With  this  scheme,  stall  margin  would  be  added 
only  for  the  duration  needed  and  amount  required  to  maintain 
suitable  engine  performance  over  the  particular  transient 
encountered . 

Active  stall  control  offers  the  potential  for  significant 
reductions  in  current  stall  margin  requirements.  Consequently, 
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sxibstantial  increases  in  fan  and  compressor  performance  and 
reductions  in  component  weight  would  be  expected.  This 
translates  into  increased  engine  thrust-to-weight  and  reduced 
specific  fuel  consximption.  Also,  improvements  in  engine 
operatbility  and  duraibility  would  be  likely  as  current  control 
schemes  could  be  designed  to  be  less  restrictive  (i.e./  greater 
throttle  movement,  fewer  afterburner  restrictions) ,  and  fewer 
compressor  stalls  encountered. 

If  any  or  all  of  the  benefits  described  above  are  to  be 
realized,  the  key  is  the  existence  and  detection  of  an  accurate 
and  robust  (able  to  be  detected  under  mild  and  severe  operating 
conditions)  stall  precursor  signal.  In  the  1960's,  limitations 
of  hydromechanical  controls  prevented  any  practical  use  of  a 
stall  precursor  signal  even  if  one  could  have  been  shown  to 
exist.  Stall  control  schemes  have  generally  relied  on  detecting 
the  presence  of  stall  (rather  than  a  stall  precursor) ,  then 
controlling  some  variable  geometry  feature  of  the  compressor  to 
lower  the  operating  point.  Ludwig  and  Nenni  [1]  demonstrated 
such  a  scheme  on  a  turbojet  engine  and  met  with  only  limited 
success,  largely  due  to  insufficient  warning  time  prior  to  the 
instability.  The  lack  of  adequate  stall  warning  time 
necessitated  the  use  of  large  control  forces  to  move  the 
compressor  operating  point. 

Recent  investigations  by  McDougall,  et  al.  [2]  and  Gamier, 
et  al.  [3]  have  identified  small  amplitude  rotating  waves  (modal 
waves)  as  a  potential  stall  precursor.  In  cases  where  the  pre- 
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stall  waves  have  been  shown  to  exist,  signals  have  been  detected 
much  sooner  than  previous  schemes,  thus  offering  potential  for  an 
active  control  system  to  take  action.  Because  of  the  substantial 
benefits  already  mentioned,  these  recent  works  have  generated  a 
lot  of  interest  in  both  the  research  and  jet  engine  manufacturing 
communities  (Aviation  Week  [4],  Hosny,  et  al  [5],  and  Paduano,  et 
al  [6]).  Although  results  to  date  have  been  promising,  several 
key  issues  must  be  resolved,  namely,  the  generality  and 
robustness  of  the  rotating  waves .  Indeed,  a  recent  investigation 
by  Day  [7]  indicated  that  modal  waves  were  not  always  present 
prior  to  stall,  and  when  present,  did  not  necessarily  show  a 
strong  coupling  to  the  stall  formation  process.  Day  [7] 
suggested  it  was  the  formation  and  siobsequent  growth  of  finite 
stall  cells  (nonlinear)  that  governed  the  stall  inception 
process.  The  works  of  McDougall,  et  al.  [2],  Gamier,  et  al. 

[3],  and  Day  [7]  are  discussed  in  more  detail  in  Chapter  2.0. 

1.2  Coapreeeor  Performance  -  Surge  and  Rotating  Stall 

As  one  of  the  major  components  of  gas  turbine  engines  (the 
others  being  the  inlet,  combustor,  turbine,  and  nozzle) ,  the 
compressor  develops  some  or  all  of  the  pressure  rise  required  by 
the  system  thermodynamic  cycle.  The  compression  system  (fan  and 
compressor)  must  operate  stablely  and  with  high  efficiency  over 
the  entire  engine  performance  envelope,  otherwise,  engine 
operability  will  be  compromised.  Overall  performance  is 
typically  displayed  on  a  compressor  map,  as  shown  in  Fig.  1.1. 
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The  compressor  operating  region  is  botinded  on  the  high  flow  side 
by  blade  passage  choking  and  on  the  low  flow  side  by  blade  and 
endwall  stalling.  Note  that  peak  pressure  rise  over  the  useful 
operating  range  of  the  compressor  occurs  right  at  the  stall  line. 
This  behavior  is  characteristic  of  most  machines  and  explains  the 
desire  to  reduce  stall  margin  requirements  (discussed  in  section 
1.3) .  Flow  breakdown  (stall)  occurs  to  the  left  of  the  stall 
line  (or  surge  line)  and  manifests  itself  as  either  of  two 
distinct  phenomena  -  surge  or  rotating  stall . 

Surge  is  a  violent  planar  dynamic  instability  characterized 
by  large  amplitude  oscillations  in  mass  flow  and  pressure 
(Fig.  1.2).  The  frequencies  of  the  oscillations  are  governed  by 
the  compressor  and  inlet  zmd  exit  (combustor)  volximes^  and  are 
usually  in  the  three  to  fifteen  hertz  range.  Because  of  the  high 
energy  levels  required  to  sustain  the  cycles,  surge  is 
encountered  at  high  compressor  speeds  (above  engine  idle) .  Most 
surges  are  self-clearing,  provided  the  initial  disturbance  has 
been  removed. 

Rotating  stall  is  a  highly  three-dimensional  localized  event 
in  which  one  or  more  packets  of  stalled  air  (stall  cells) 
propagate  circumferentially  around  the  conqpressor  at  speeds 
between  twenty  and  seventy  percent  of  rotor  speed  (Fig.  1.2)  . 

The  circumferential  extent  of  the  stall  cells  can  range  from  only 
a  few  blade  passages  to  over  180  degrees  of  the  annulus.  The 
cells  can  be  full  or  part  spem  and  may  or  may  not  cover  the  full 
axial  extent  of  the  compressor.  Rotating  stalls  can  be 


4 


classified  as  abrupt  or  progressive,  in  relation  to  the 
performance  behavior  as  viewed  on  a  compressor  map.  Figure  1.2b 
is  schematically  representative  of  an  abrupt  stall,  evidenced  by 
the  sudden  discontinuous  decrease  in  performance  once  stall 
occurs .  Progressive  stall  is  characterized  by  a  smooth, 
continuous  decrease  in  performance,  typical  of  a  part-span  stall 
which  grows  slowly  in  circumferential  and  radial  extent  (Greitzer 
[12]) ,  Both  types  of  rotating  stall  (abrupt  and  progressive) 
were  observed  in  the  current  set  of  experiments. 

Rotating  stall  is  characterized  by  quasi-steady  (globally 
stable  in  a  mass-averaged  sense)  performance  at  a  reduced  flow 
rate  and  pressure  rise.  It  is  generally  accompanied  by  recovery 
hysteresis;  that  is,  once  rotating  stall  is  fully  developed, 
recovery  to  normal  operation  requires  a  significant  increase  in 
exit  throttle  area  (applicable  in  test  cells;  note  in  engines, 
bleed  flow,  variable  vanes  and/or  fuel  control  must  be  used) 
beyond  that  at  which  stall  inception  occurred  (Fig.  1.2b) . 
Extensive  hysteresis  can  prevent  recovery  altogether,  thus 
leading  to  a  "hung"  stall.  This  condition  has  become  known  as 
stagnation  or  non-recoverable  stall. 

In  aircraft  gas  turbine  engines,  surge  is  the  predominant 
compressor  instability  encountered.  However,  continuous  surge 
cycles  (non-clearing)  result  in  decreasing  rotor  speed  which  can 
cause  the  instability  to  degenerate  into  a  rotating  stall  and 
possible  non-recoverable  condition  (Benrey  [8]).  In  either  case, 
consequences  are  damaging  and  can  be  severe.  Both  surge  and 
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rotat.ing  stall  lead  to  large  vibratory  stresses  in  the  compressor 
blading  aind  possible  turbine  inlet  overtenperature . 

Additionally,  the  reduced  performance  levels  associated  with 
rotating  stall  can  result  in  an  engine  cycle  that  is  not  self- 
sustaining  . 

It  is  worth  mentioning  that  most  engine  conqpressor  stalls 
are  caused  from  a  highly  transient  event  (hard  afterburner  light, 
severe  inlet  distortion  from  a  sudden  maneuver,  hot  gas  ingestion 
from  an  armament  firing,  etc.).  In  a  compressor  rig  test,  stall 
is  typically  induced  by  increasing  the  backpressure  via  a 
downstream  throttle  (as  was  done  in  the  current  investigation) . 
Thus,  from  a  practical  standpoint,  it  is  useful  to  study  the 
effects  of  different  throttle  closure  rates  on  the  stalling 
behavior  of  compressor  rigs  (faster  rates  more  closely  associated 
with  actual  engine  stalls) .  In  the  present  experiments,  only  one 
throttle  rate  was  available  (discussed  in  detail  in  section  4.1). 
The  effects  of  throttle  closure  rate  are  discussed  further  in 
section  2.4. 


1.3  Stability  Audit  -  Stall  Margin 

The  need  to  avoid  the  severe  conseq[uences  of  rotating  stall 
and  surge  inqposes  a  significant  dilemma  on  the  engine  designer; 
the  engine  must  optimally  satisfy  the  conflicting  requirements  of 
high  thrust,  low  fuel  consumption,  low  weight,  long  life,  low 
cost,  and  adequate  engine  stability.  The  stability  audit  (or 
stackup)  is  a  procedure  currently  used  to  establish  required 
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surge  margins .  Figure  1 . 3  shows  the  stability  audit  items  or 
destabilizing  influences  that  must  be  included  in  determining  the 
design  stall  margin  of  a  high  pressure  compressor.  Many  of  the 
factors  are  not  usually  present  (or  not  present  all  at  the  same 
time) ;  consequently,  the  compressor  operates  with  more  stability 
margin  than  is  required.  Typical  design  values  of  stall  margin 
range  from  15-25%  (Wisler  [9]),  defined  as 

SM-  (1.1) 


at  a  constant  flow. 

A  signal  that  consistently  precedes  a  stall  event  (rotating 
stall  or  surge)  with  enough  warning  time  to  take  action  to  avoid 
or  suppress  the  instability  could  significantly  reduce  the  surge 
margin  requirements  of  fan  and  compressor  components.  Ideally, 
if  a  robust  stall  precursor  did  exist,  active  controls  would  be 
used  to  interfere  with  the  unsteady  aerodynamic  damping  of  the 
compressor  in  a  manner  that  moved  the  point  of  instability  onset 
to  lower  flows.  This  idea  of  an  actively  stabilized  stall  line 
was  first  proposed  in  the  open  literature  by  Epstein,  et  al. 

[10],  and  is  illustrated  in  Fig.  1.4.  In  Fig.  1.4,  a  compressor 
is  allowed  to  accelerate  at  a  maximum  rate  (transient  operating 
line)  and  actually  operate  to  the  left  of  the  steady-state  stall 
boundary  (net  stall  line) .  Once  the  acceleration  is  complete, 
the  additional  margin  is  not  required  and  the  stall  line  is 
returned  to  the  base  stall  line.  Thus,  the  compressor  operating 
line  is  raised  and  adequate  stall  margin  provided  only  for  the 
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duration  and  amount  needed  to  maintain  stable  performance  over 
the  transient  event.  The  work  of  Paduano,  et  al.  [6]  and  Day 
[11]  are  two  recent  exaunples  of  the  use  of  active  compressor 
stabilization . 

1.4  Intent  of  Current  Work 

The  primary  objective  of  the  current  investigation  is  to 
characterize  the  stall  inception  behavior  in  two  high-speed, 
advanced  design  compressors .  Results  are  intended  to  provide  a 
sound  physical  understanding  of  the  stall  development  process  in 
machines  that  are  typical  of  those  used  in  current  zmd  future 
aircraft  engine  applications .  It  is  emphasized  that  the  majority 
of  reported  stall  investigations  performed  to  date  have  been 
carried  out  on  low-speed,  low  pressure  rise  machines 
(compressibility  effects  are  assumed  negligible) .  Questions 
regarding  the  appliceUsility  of  these  results  to  high-speed 
(compressible)  machines  must  be  answered.  Thus,  the  current 
stall  investigation  work  attempts  to  provide  information  that 
will  be  useful  in  assessing  this  applicability. 

The  two  single-stage  rotor  designs  differ  primarily  in  their 
amount  of  leading  edge  sweep  (one  has  a  conventional  straight 
leading  edge,  the  other  is  swept  back) .  Both  were  tested  at 
several  speeds  under  various  clean  inlet  (undistorted)  operating 
conditions.  In  addition  to  numerous  instrumentation  to  quantify 
the  compressor  steady  performance,  the  rig  case  was  instnunented 
with  eight  high-response  pressure  transducers  equally  spaced 
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around  the  annulus  for  stall  development  detection.  The  data 
were  analog  recorded,  then  post-processed  as  the  intent  of  the 
current  work  is  to  characterize  the  stall  inception  behavior,  and 
not  to  demonstrate  a  real-time  control  strategy. 

A  critical  question  that  must  be  answered  is  which  process, 
modal  waves  or  finite  cells  (last  paragraph,  section  1.1), 
dominates  the  formation  of  fully  developed  stall  (rotating  stall 
or  surge)  .  Resolution  of  this  issue  is  key  if  inqprovements  in 
the  modelling  of  compressor  stability  are  to  be  realized  and  the 
most  beneficial  approaches  to  active  control  defined.  The 
current  work  focuses  on  this  issue,  and  in  particular  addresses 
the  following  questions : 

(1)  What  is  the  relationship  between  modal  waves  and  finite 
cells,  and  is  the  relationship  consistent  or  does  it  vary  with 
compressor  design,  operating  condition,  etc? 

(2)  What  is  the  effect  of  compressor  operating  condition  on 
the  pre-stall  indication? 

(3)  Where  precisely  on  the  compressor  pressure  rise 
characteristic  do  the  first  signs  of  pre-stall  behavior  occur? 
Does  the  shape  of  the  characteristic  affect  the  pre-stall 
behavior? 

(4)  How  do  the  features  of  the  stall  inception  process  in 
high-speed  machines  compare  to  those  in  low-speed  machines? 

(5)  Are  there  potential  stall  warning  signs  that  might  be 
peculiar  to  high-speed  compressors? 

Additionally,  a  whole  series  of  questions  regarding  the 
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robustness  of  the  pre-stall  indications  (influence  of  distorted 
flow  fieldSf  various  throttle  rates,  etc.)  must  be  addressed,  but 
these  are  not  exeunined  in  this  work. 
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2 . 0  BACKGROUND 


2 . 1  Preliminary  Remarks 

As  noted  by  Greitzer  [12],  the  subject  of  axial  con^ressor 
stall  (rotating  stall  etnd  surge)  can  be  separated  conceptually 
into  two  main  areas: 

(1)  Stall  inception  -  This  area  deals  with  the  examination 
of  the  basic  fluid  mechanics  associated  with  the  onset  of  stall. 
Important  aspects  include  the  influence  of  inlet  distortion  and 
rapid  throttle  transients,  the  effects  of  downstream  system 
components,  and  stability  enhancements  such  as  rotor  casing 
treatments  and  active  controls. 

(2)  Post-stall  behavior  -  In  this  area,  behavior  subsequent 
to  the  onset  of  the  initial  instability  is  of  prime  concern. 
Topics  of  interest  include  large  an^litude  surge  oscillations, 
transient  performance  into  and  out  of  rotating  stall,  structure 
of  the  rotating  stall  flow  field  and  in-stall  operation,  and 
recovery  hysteresis  extent. 

Because  of  the  nature  of  the  current  work  this  literature  survey 
will  focus  on  the  first  area,  stall  inception.  In  particular, 
those  theories  and  investigations  that  relate  to  the  two 
currently  popular  explanations  of  stall  development,  modal  waves 
and  finite  cells,  will  be  enq>hasized.  Excellent  summaries  of 
previous  work  in  the  area  of  post- stall  performance  are  provided 
by  Davis  [13],  Copenhaver  [14],  and  Gorrell  [15]. 

As  was  mentioned  in  section  1.4,  the  majority  of  stall 
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investigations  performed  to  date  (both  analytical  and 
experimental)  have  been  carried  out  on  low-speed,  high  aspect 
ratio  machines  with  limited  niamber  of  stages  (low  presure  rise)  . 
The  theories  that  have  developed  from  these  investigations  have 
been  able  to  add  considerable  insight  into  the  fluid  mechanics  of 
stability  onset  and  stall  behavior.  However,  they  have  limited 
use  in  regards  to  the  high-speed,  high  pressure  ratio  aircraft 
engine  applications  which  are  of  practical  interest.  The 
limitations  of  the  theories  are  imposed  by  the  assumptions  that 
are  inherent  in  their  development,  one  of  the  most  crucial  being 
incompressible  treatment  of  the  flow.  The  addition  of 
compressibility  adds  much  complexity  to  the  theories,  but  is 
essential  if  practical  use  is  to  come  from  the  investigations. 

The  work  that  has  been  performed  on  high-speed  compressors  has 
largely  been  involved  with  post -stall  and  recovery  behavior. 
Recent  examples  (1985-1991)  include  Davis  [13],  Copenhaver  [14], 
Gorrell  [15],  Small  and  Lewis  [16],  Hosny  and  Steenken  [17], 

Boyer  and  O'Brien  [18],  [19],  and  Bonnaure  [20].  Of  these,  only 
Bonnaure's  work  [20]  focused  on  the  stall  inception  process 
through  the  development  of  a  two-dimensional  model  of  high-speed 
compressor  stability.  The  present  research  focuses  on  an 
experimental  investigation  of  the  stall  development  process  in 
high-speed,  advanced  design  con^ressors . 

2 . 2  Karly  Korks 

For  the  purpose  of  this  literature  review,  work  dealing  with 
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stall  inception  performed  prior  to  1980  is  classified  as  early. 

An  excellent  review  of  these  earlier  analyses  is  provided  by 
Greitzer  [12] . 

2.2.1  Basic  Mechanism  of  Rotating  Stall 

A  simple  and  still  widely  used  explanation  of  the  mechanism 
of  rotating  stall  was  initially  proposed  by  Emmons,  et  al  [21] 
and  recently  applied  by  Day  [7] .  Consider  a  row  of  blades 
operating  on  the  verge  of  stall  (large  angle  of  attack)  as  in 
Fig.  2.1.  If  some  flow  disturbance  causes  one  (or  several 
adjacent)  of  the  blades  to  momentarily  stall,  then  the  flow 
through  that  passage  will  be  restricted  by  the  separated  boundary 
layer  on  the  suction  side  of  the  stalled  blade  (blade  A  of  Fig. 
2.1) .  Consequently,  the  incoming  flow  will  be  diverted  to 
surrounding  passages,  increasing  the  incidence  on  blade  B  and 
decreasing  the  incidence  on  blade  C.  Thus,  the  stall  will 
propagate  from  blade  to  blade  in  a  direction  opposite  of  rotation 
relative  to  the  blades,  but  in  the  direction  of  rotation  in  em 
absolute  reference  frame.  If  the  initial  flow  disturbance  is  not 
relieved  or  is  worsened,  the  stall  can  grow  into  a  fully 
developed  cell  covering  more  than  half  of  the  compressor  annulus. 
Experimental  evidence  has  shown  rotating  stall  cells  to  propagate 
anywhere  between  twenty  and  seventy  percent  of  the  rotational 
speed  of  the  blades,  depending  on  the  size  and  extent  of  the 
cell . 

The  Emmons  model  (Fig.  2.1)  thus  suggests  a  possible  cause 
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of  rotating  stall  inception;  namely,  an  instability  associated 
with  the  stalling  of  a  blade  passage <s) .  Day  [7]  recently  used 
this  flow  model  to  offer  an  explanation  of  stall  inception  based 
on  disturbances  of  a  short  length  scale  (a  few  blade  passages) . 
Day's  work  is  discussed  in  more  detail  in  section  2.3.2. 

2.2.2  Linearized  Theories  of  Stall  Inception 

The  basis  for  the  modal  wave  explanation  of  stall 
development  stems  from  hydrodynamic  stability  analysis,  which 
deals  with  small  auaplitude  disturbances  and  thus  linearized  flow 
equations .  The  stability  issue  involves  examination  of  the 
growth  or  decay  of  the  disturbances.  Rotating  stall  and  surge 
are  viewed  as  limit  cycles  whose  final  strength  is  governed  by 
non-linear  effects. 

The  linearized  stability  analyses  for  rotating  stall  onset 
share  many  common  features  in  their  treatment  of  the  compressor 
flow  field.  Typically,  the  models  are  two-dimensional  and 
incompressible,  and  describe  the  evolution  of  a  small  amplitude 
perturbation  superimposed  on  a  mean  compressor  operating 
condition.  The  wavelength  of  the  perturbations  is  taken  to  be 
much  larger  than  the  blade  pitch  so  that  the  blade  rows  are 
represented  as  actuator  or  semi-actuator  disks^.  Conditions 
applied  across  the  blade  row  to  link  the  flow  quantities  upstream 

^Defined  as  a  model  of  a  blade  row  as  a  plane  across  which 
mass  flow  is  continuous,  but  other  properties  can  be 
discontinuous.  The  semi-actuator  model  attempts  to  represent 
con^ressor  inertial  effects,  usually  through  a  first-order  lag 
equation  (Greitzer  [12]). 
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and  downstream  are  mass  conservation,  an  inlet/exit  flow  angle 
relation,  and  a  momentum  relation.  Differences  between  the 
models  are  usually  memifested  in  the  form  of  the  momentum 
relation  used  by  each.  These  forms  have  included  static  pressure 
rise  (Stenning  [22]),  total  pressure  loss  (Fabri  [23]),  and 
vorticity  conservation  with  friction  losses  (Nenni  and  Ludwig 
[24]),  all  described  as  functions  of  inlet  flow  angle.  The 
linearized  equations  of  motion  are  solved  to  determine  the 
eigenvalues  which  define  the  stability  of  the  flow  field,  or  to 
examine  the  growth  of  the  initially  prescribed  small  perturbation 
(until  the  disturbance  is  too  large  that  the  linearized  treatment 
is  no  longer  valid) . 

All  of  the  theories  require  blade  row  or  stage  performance 
data  as  input,  which  typically  must  be  determined  by 
experiment^.  Additionally,  a  stability  criterion  to  define  the 
onset  of  stall  development  must  be  specicifed.  This  is  usually 
defined  to  be  the  neutral  stability  criterion,  taken  as  the  zero 
slope  point  on  the  compressor  total-to-static  pressure 
characteristic  (Fig.  2.2).  As  pointed  out  by  Greitzer  [12],  this 
zero  slope  condition  provides  a  rough  "rule  of  thumb, "  but  is  not 
applicable  to  all  machines. 

Despite  their  limitations,  the  linearized  theories  have  been 
useful  in  improving  the  understanding  of  the  axisymmetric  flow 
breakdown  process  in  axial-flow  compressors.  Additionally,  they 

^The  prediction  of  compressor  performance  characteristics  at 
off-design  conditions  is  difficult,  and  at  present,  cannot  be 
done  with  the  required  accuracy. 
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have  been  used  to  study  the  stage  (and/or  blade  row)  interaction 
and  overall  compressor  system  effects  on  staJbility.  These 
interactions  alter  the  stage  stalling  behavior  from  that  of  the 
ssune  stage  operating  in  isolation.  These  areas  continue  to  be  of 
prime  interest  (Longley  and  Hynes  [25] ) . 

2 . 3  Recent  Works 

As  noted  in  the  Introduction  of  this  report^  two 
explanations  of  the  compressor  stall  inception  process  are 
currently  popular;  small  an^litude  rotating  waves  (modal  waves) 
emd  short  length  scale  finite  disturbances  (finite  cells) .  This 
section  of  the  literature  survey  will  examine  these  two  flow 
models  in  more  detail.  Additionally,  recent  experimental 
applications  of  both  models  in  the  area  of  active  suppression  of 
rotating  stall  will  be  briefly  discussed. 

2.3.1  Modal  Waves  as  a  Stall  Inception  Indication 

The  wave  model  of  flow  field  dynamics  evolved  partly  from 
the  early  theories  described  above  and  more  directly  from  the 
theory  of  rotating  stall  transients  proposed  by  Moore  [26]  and 
Moore  and  Greitzer  [27] .  In  their  rotating  stall  2uialysis,  Moore 
and  Greitzer  represent  the  upstream  flow  coefficient,  ^  ,  as 
consisting  of  a  steady  mean  flow  consonant,  ^  ,  and  a  non- 
axi symmetric  disturbance,  : 

(2.1) 

Gamier  [28]  took  the  equation  which  models  the  evolution  of  the 
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non-axisyinmetric  component  and  linearized  it  for  the  case  of 
small  an^litude  disturbances.  The  resulting  governing  equation 
for  the  non-axisymmetric  disturbance  was  shown  to  be  of  the  form: 

2  (2.2) 

k*0 

where:  k  «  a  Fourier  mode  (k>0  corresponds  to  surge-like 

disturbances ) 

(0  »  wave  frequency  (according  to  theory.  Gamier 
[28],  should  correspond  to  rotating  stall 
frequency) 

a  •  growth  rate  or  damping  of  the  wave 
6  «  circximferential  position 
t  «  time 

In  essence,  the  modal  wave  view  treats  the  compressor  system 
as  a  harmonic  oscillator  with  stall  onset  associated  with  the 
growth  of  long  wavelength  (long  compared  to  blade  pitch) 
circTomferential  disturbances .  Surge  and  rotating  stall  are 
viewed  as  eigenmodes  of  the  system  (Fig.  2.3),  with  surge  being 
of  zero  order  and  rotating  stall-like  disturbances  being  higher 
order.  For  example,  a  velocity  perturbation  with  a  wavelength 
equal  to  the  circumferential  length  of  the  compressor  annulus 
would  be  considered  a  mode  of  order  one.  The  growth  or  decay  of 
the  instabilities  depend  upon  the  stability  criteria  related  to 
the  operating  point  of  the  conqsressor.  If  the  compressor  were 
throttled  towards  the  stall  point,  the  stall  cell  would  grow 
smoothly  out  of  the  small  amplitude  modal  waves. 
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Recent  experimental  investigations  by  McDougall,  et  al.  [2] 
and  Gamier,  et  al.  [3]  have  provided  support  to  the  wave  model 
of  stall  inception.  McDougall,  et  al.  [2]  were  the  first  to 
demonstrate  the  rotating  disturbances  on  a  single-stage,  low- 
speed  axial  compressor.  Two  methods  were  used  for  detection: 

(1)  cross-correlation  of  the  unsteady  part  of  signals  from  two 
hot  wires  set  90  degrees  apart  0 . 6  mean  radii  upstream  of  the 
rotor,  and  (2)  spatial  and  temporal  Fourier  euialysis  of  signals 
simultaneously  logged  from  six  hot  wires  equally  spaced  around 
the  circumference,  again  0.6  radii  upstream.  Resolution  of  the 
circumferential  harmonics  showed  the  first  order  Fourier  mode  to 
be  the  dominant  one,  rotating  in  synqpathy  with  the  variation  in 
endwall  blockage  at  rotor  exit.  The  small  amplitude  upstream 
perturbations  (about  0.5%  of  the  mean  axial  velocity  near  the 
stall  point)  were  shown  to  rotate  around  the  annulus  at  about  49% 
of  rotor  speed  and  tremsition  into  fully  developed  rotating  stall 
with  no  discontinuity  in  phase  or  frecpiency.  The  waves  could  be 
detected  as  much  as  two  seconds  prior  to  stall,  rising  and 
falling  in  amplitude,  but  were  of  significant  amplitude  only  at 
conditions  very  close  to  stall  (one  rotational  speed  was  examined 
and  the  stall  was  initiated  in  a  quasi-steady  manner) . 

McDougall,  et  al.  [2]  concluded  that  the  stalling  process  was 
closely  related  to  large  length  scale  disturbances  (lower  order 
modes)  and  not  to  events  occurring  in  individual  blade  passages. 

Gamier,  et  al.  [3]  detected  rotating  waves  in  two  low-speed 
compressors  (a  single-stage  and  a  three-stage)  and  one  high- 
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speed,  three-stage  compressor.  The  compressors  were  operated  in 
uniform  inlet  flow,  distorted  inlet  flow,  and  subject  to  various 
throttle  rates  at  a  single  rotational  speed.  Discrete  Fourier 
analysis  of  the  signals  was  used  to  extract  the  modal  information 
(hot  wires  were  used  for  the  low-speed  rigs,  high  response  static 
pressure  transducers  for  the  high-speed  machine) .  Sets  of  high- 
response  instrvunentation  were  mounted  at  several  axial  stations 
to  exzunine  sensor  placement  influence.  Again,  the  waves  were 
initially  of  small  amplitude  and  evolved  into  rotating  stall 
without  sharp  changes  in  phase  or  frequency.  The  signals  were 
clearest  when  measured  by  the  sensors  nearest  to  the  stage  which 
initiated  the  stall.  Figure  2.4  presents  a  szunple  of  Gamier' s, 
et  al.  low-speed  results  [3]  using  the  Fourier  coefficients.  The 
amplitude  of  the  first  Fourier  harmonic  is  shown  to  be  non- zero 
for  2a3out  ninety  rotor  revolutions  prior  to  stall.  Additionally, 
the  phase  speed  of  the  first  harmonic  is  essentially  constant 
over  the  same  period,  and  approximately  equal  to  the  rotating 
stall  speed  (35%  rotor  speed  before  stall,  38%  speed  in  fully 
developed  stall).  Note  as  with  McDougall's,  et  al.  work  [2],  the 
first  order  mode  was  dominant  in  the  low-speed  results.  Results 
from  Gamier' s,  et  al.  high-speed  investigations  [3]  were 
qualitatively  similar  to  the  low-speed  results,  but  also  revealed 
significant  influence  by  the  second  harmonic.  No  explanations 
have  been  offered  for  the  dominance  of  one  mode  over  another. 

Gamier,  et  al.  [3]  concluded  that  the  modal  waves  and  fully 
developed  rotating  stall  are  two  stages  of  the  same  phenomenon. 
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He  also  noted  that  simple  Fourier  techniques  were  inadequate  for 
determination  of  the  waves  in  distorted  flow  fields,  and 
suggested  a  signal  processing  method  based  on  the  true  system 
eigenmodes . 

2.3.2  Finite  Cells  as  a  Stall  Inception  Indication 

Day' s  recent  experimental  work  with  two  low- speed 
compressors  [7,  11]  has  resulted  in  an  alternative  view  of  the 
stall  development  process;  namely,  the  importance  of  small  length 
scale  disturbances  (finite  cells) .  Using  similar  instrumentation 
and  data  processing  as  McDougall,  et  al.  [2]  and  Gamier,  et  al. 
[3] ,  Day  [7]  showed  that  modal  waves  did  not  always  precede 
stall,  and  when  present  did  not  necessarily  show  a  strong 
coupling  to  the  stall  development  process.  Day's  work  showed  the 
stall  inception  process  to  be  governed  by  the  formation  and 
sxibsequent  growth  of  finite  cells  of  length  scales  he  classified 
as  small  and  large  with  respect  to  blade  pitch.  Increased 
coupling  of  the  modal  waves  and  finite  cells  was  more  likely  when 
the  length  scales  of  the  waves  and  stall  cell  were  comparable 
(i.e.,  when  the  stall  cell  covered  a  large  region  of  the 
annulus) .  Tip  clearance  appeared  to  be  important  in  determining 
which  stall  inception  model  applied,  with  increased  clearance 
more  likely  to  reveal  modal  waves  prior  to  the  development  of  a 
longer  length  scale  finite  cell  (Fig.  2.5) .  Tighter  tip 
clearance  favored  the  formation  of  a  smaller  length  scale  cell 
with  no  evidence  of  modal  waves  beforehand  (Fig.  2.6) .  Day 
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concluded  that  modal  perturbations  and  the  formation  of  finite 
stall  cells  were  not  necessarily  related;  either  disturbance 
might  be  the  first  to  appear. 

2.3.3  Active  Control  experiments 

Although  not  directly  related  to  the  scope  of  the  present 
work^  recent  experiments  by  Paduano^  et  al.  [6]  and  Day  [11] 
using  active  controls  to  suppress  rotating  stall  and  delay  its 
onset  demonstrate  the  potential  benefits  of  identifying  a 
reliad>le  stall  precursor  signal.  Paduano,  et  al.  [6]  applied  the 
modal  wave  flow  model  by  "wiggling"  inlet  guide  vanes  to  danqp  the 
perturbations  in  the  same  single-stage  compressor  used  by 
Gamier,  et  al.  [3].  This  resulted  in  a  20%  gain  in  compressor 
mass  flow  range  when  the  first  two  Fourier  modes  were  controlled. 

Day  [11]  used  fast -acting  air  injection  valves  both  globally 
(the  entire  aimulus)  to  damp  out  modal  waves  (when  present)  and 
locally  (a  portion  of  the  annulus)  to  remove  emerging  finite 
stall  cells .  Results  indicated  a  4%  increase  in  flow  range  when 
the  waves  were  damped  and  a  6%  increase  in  range  when  the  finite 
cells  were  removed  (the  valves  could  not  be  operated  in  both 
global  and  local  modes  simultaneously).  Both  Paduano's,  et  al. 

[6]  and  Day's  [11]  active  control  schemes  required  minimal  power 
to  inclement  conqpared  to  the  power  of  the  compressor  (a  major 
difference  from  previous  stall  control  schemes) ,  and  demonstrated 
the  feasibility  of  actively  suppressing  compressor  instabilities. 


21 


2.4  Potential  Factors  Affecting  Stall  Precursor 

It  is  appropriate  to  provide  some  discussion  regarding  the 
influence  (or  potential  influence)  of  a  variety  of  factors  on 
stall  precursor  signals  (or  the  aibility  to  detect  such  signals)  . 
These  factors  include  compressor  environment  (inlet  distortion) , 
operating  conditions  (rotational  speedy  mass  flow  transients) , 
sensor  location,  speedline  shape  and  factors  which  influence  that 
shape  (i.e.,  tip  clearance),  and  compressibility  effects  (shock 
waves,  stage  mismatching,  etc.).  Investigating  the  influence  of 
these  factors  is  important  for  developing  a  so\and  physical 
understanding  of  the  stall  inception  process  so  that  the  factors 
(or  their  effects)  can  be  included  in  advanced  modelling  of 
compressor  stedsility. 

Several  of  these  factors  have  been  documented  in  previous 
investigations.  Gamier  [28]  showed  that  inlet  distortion 
(spatially  non-uniform  inlet  total  presure)  hindered  detection  of 
the  pre-stall  waves  (using  Fourier  techniques) .  He  further 
suggested  an  analysis  technique  based  on  the  true  system 
eigenmodes  (independent  of  wave  shape) .  The  effect  of  sensor 
location  was  also  studied  by  Gamier  [28]  in  his  single-stage, 
low-speed  and  three-stage,  high-speed  compressor  investigations. 
The  low-speed  results  indicated  no  major  differences  in  terms  of 
when  wave  propagation  could  be  detected  versus  axial  position  of 
sensors.  The  high-speed  results  showed  much  stronger  dependence 
on  axial  sensor  location,  with  maximvim  stall  warning  time 
obtained  from  the  measurements  (leading  edge  of  the  first  stator) 


closest  to  the  stage  which  initiated  the  stall.  Additionally/ 
Gamier  [28]  examined  the  effect  of  different  throttle  closure 
rates  (mass  flow  transients) /  and  found  the  stall  warning  time  to 
be  roughly  inversely  proportional  to  the  time  rate  of  change  of 
mass  flow  at  stall.  It  is  importmt  to  note  that  the  throttle 
closure  rates  used  by  Gamier  [28]  (maximiim  d({>/dt~-0 . 00073  1/s) 
all  resulted  in  quasi-steady  compressor  stalling  (i.e.,  the 
compressor  steady  stall  line  was  valid) .  Rapid  closure  rates 
would  likely  result  in  transient  compressor  stalling  (steady 
stall  line  no  longer  valid) ,  discussed  further  in  section  4.1. 

The  importance  of  speedline  shape  on  compressor  stalling 
behavior  has  been  well  documented  in  several  experimental  and 
analytical  investigations  (i.e.,  Greitzer  [31]/  Bonnaure  [20]). 

In  fact/  all  of  the  compressor  post-stall  and/or  stability  models 
rely  on  accurate  compressor  characteristics  (speedlines)  to 
represent  the  physical  behavior.  As  discussed  in  section  2.2.2/ 
the  zero  slope  condition  is  typically  taken  as  the  onset  of 
compressor  unstable  performance.  For  advanced  stall  warning 
time/  a  characteristic  with  a  local  positive  slope  at  stall  might 
be  desirable.  A  speedline  that  is  bent  over  (positively  sloped) 
is  identified  with  thickening  boundary  layers/  larger  wakes /  and 
in  the  extreme/  endwall  or  tip  stalling.  These  are  all  features 
of  impending  flow  breakdown  -  stall.  Note  that  tighter  tip 
clearances  would  tend  to  reduce  the  positively  sloped  region  of 
the  characteristic  (less  losses)  and  result  in  potentially  leas 
stall  warning.  Indeed/  Day's  results  [7]  showed  that  increased 
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tip  cleareuice  appeared  to  favor  formation  of  modal  waves  which 
evolved  smoothly  into  rotating  stall  (see  Fig.  2.5). 

The  influence  of  compressibility  effects  on  con^ressor  pre¬ 
stall  behavior  has  not  received  as  much  attention  as  the  before- 
mentioned  factors  (largely  due  to  the  lack  of  high-speed 
experimental  results) .  In  his  recent  analytical  investigation, 
Bonnaure  [20]  found  the  influence  of  compressible  unsteady 
perturbations  (local  mass  flow  perturbations  which  vary  along  the 
compressor  length;  in  low-speed  case,  these  are  the  same 
everywhere)  to  be  key  in  determining  compressor  stability. 

Indeed,  Bonnaure' s  work  [20]  appeared  to  confirm  the  high-speed 
results  of  Gamier  [28]  which  showed  strong  dependence  on  axial 
sensor  location  (preceding  paragraph) . 

The  effect  of  shock  waves  on  the  stalling  process  of 
transonic  designs  (as  are  those  in  the  current  investigation)  is 
not,  in  general,  well  luiderstood.  For  exanqple,  much  work  is 
currently  ongoing  in  the  facility  where  the  current  tests  were 
conducted  (see  section  3.1)  to  characterize  the  rotor  tip  shock 
structure  of  several  designs  at  near-stall  operating  conditions. 
Although  the  current  set  of  experiments  did  not  quantify  the 
rotor  shock  strength,  the  one  significant  design  difference 
(backwards  sweep  versus  straight  leading  edge)  allowed  some 
inferences  to  be  made  regarding  the  potential  influence  of  shock 
waves  on  stall  precursor  signals.  A  discussion  of  this  influence 
is  provided  in  section  5.6. 
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3.0  TS8T  APPARATUS 


3 . 1  Facility  Description 

The  facility  used  for  all  the  testing  was  the  Aero 
Propulsion  and  Power  Directorate's  2000  horsepower  (hp) 

Compressor  Aerodynamic  Research  Laboratory  (CARL)  located  at 
Wright -Patterson  AFB,  Ohio.  An  artist's  rendition  of  the 
facility  is  illustrated  in  Fig.  3.1.  Table  3.1  summarizes  the 
test  cell  operating  parameters.  The  CARL  is  a  closed  loop 
facility  shown  schematically  in  Fig.  3.2. 

Air  passes  through  the  30-inch  diameter  inlet  duct  to  a 
Universal  venturi  tube  located  six  pipe  diameters  downstream  of 
the  return  tube  elbow.  The  air  is  turned  90-deg  with  the  aid  of 
turning  vanes  and  passes  through  a  ttibe  bxindle  into  a  48-inch 
diameter  settling  chamber.  The  screens  installed  both  upstream 
and  downstream  of  the  turning  vanes  are  used  to  prevent  feedback 
related  to  flow  separation  on  the  vanes  from  reaching  the 
venturi.  From  the  settling  chamber,  air  enters  the  compressor 
through  a  direct -coupled  bellmouth.  The  flow  leaving  the 
con^ressor  is  deflected  radially  outward  to  a  peripheral 
throttle.  The  throttle  consists  of  one  stationary  and  one 
rotating  cylindrical  ring,  each  with  16  circumferentially 
distributed  holes.  The  throttle  is  designed  to  vary  continuously 
from  fully  open  to  fully  closed,  and  throttling  takes  place  at  a 
dieuneter  of  eUsout  47  inches .  A  surge  valve  is  opened  to  bypass 
the  throttle  when  recovering  from  conqpressor  stall .  Downstream 
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of  the  throttle,  the  air  enters  a  collector  and  then  passes 
through  a  24-inch  diameter  duct  to  a  heat  exchanger  and  filter. 
The  air  is  filtered  to  remove  5-micron  particles  and  then 
returned  to  the  facility  through  the  30-inch  diameter  inlet  duct. 

3.2  Test  Vehicle  and  Designs  Tested 

One  of  the  primary  missions  of  the  CARL  is  to  conduct 
Paraunetric  Blade  Studies  (PBS)  to  investigate  the  effects  of 
specific  rotor  blade  design  parameters  (maximvim  thickness 
location,  surface  angles,  leading  edge  sweep,  etc.)  on  the 
performance  of  one  compressor  configuration.  To  date,  ten  rotor 
designs  euid  one  baseline  have  been  tested  under  the  PBS  series . 
For  the  current  research,  two  of  these  rotors  were  used.  This 
section  provides  brief  descriptions  of  the  rig  configuration  and 
rotors  tested  specifically  for  the  present  work.  More  general 
information  regarding  the  PBS  tests  can  be  found  in  a  report  by 
Law  and  Puterbaugh  [29] . 

3.2.1  Fan  Rig  Description 

A  cross-section  of  the  rig  is  shown  in  Fig.  3.3.  The  rotor 
tip  diameter  is  constant  at  17  inches.  The  configuration  used  no 
inlet  guide  vanes.  The  rotor  shaft  is  movmted  on  an  cil-dasiped 
roller  bearing  at  the  forward  location  and  a  ball  bearing  at  the 
aft  location.  The  rotors  were  machined  from  single  forgings’  of 

’These  single  castings  are  commonly  referred  to  as  blisks 
for  bladed  disks,  or  IBRs  for  integrally  bladed  rotors. 
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6A1-4V  titani\jia.  The  stator  was  also  fabricated  as  an  integral 
ring  machined  from  AMS  5616.  The  backward  swept  stator  design 
used  a  controlled  diffusion  airfoil  (CDA)  approach  for 
controlling  suction  surface  diffusion  to  prevent  boxindary  layer 
separation.  Stator  design  parameters  are  summarized  in  Table 
3.2. 

3.2.2  Rotor  Designs  Tested 

Rotors  4  and  6  of  the  PBS  series  were  investigated  for  the 
current  stall  precursor  research.  Table  3.2  provides  some  of  the 
design  features  of  the  rotors.  Both  are  considered  low-aspect 
ratio,  high-throughflow,  transonic  designs.  At  design 
conditions,  rotor  tip  speed  is  1500  ft/s,  and  approximately  60% 
of  the  blade  span  sees  supersonic  relative  flow  velocities.  The 
leading  edge  of  Rotor  6  is  swept  backwards  to  reduce  shock 
losses.  Rotor  4  has  a  more  conventional  straight  leading  edge. 
Photographs  of  the  rotors  installed  in  the  rig  are  shown  in  Figs. 
3.4  and  3.5.  The  location  of  the  high-response  pressure 
transducers  (Kulites)  relative  to  the  leading  edge  tip  of  each 
rotor  is  also  indicated.  Close-ups  of  this  area  (Figs.  3.6  and 
3.7)  reveal  that  the  distance  from  the  blade  tips  to  the  Kulites 
was  the  same  for  both  rotors,  2.5  inches. 

Tip  clear2mce  as  a  function  of  rotor  speed  for  Rotors  4  and 
6  is  shown  in  Fig.  3.8.  Note  that  the  clearances  are  nearly 
identical  for  both  designs.  Hot  (rotating)  clearance  was 
measured  with  an  active,  non-touching,  spark-gap  type  clearance 
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system  at  the  rotor  leading  edge  and  mid-chord  regions  (encircled 
areas  1  and  2,  respectively,  in  Fig.  3.3)  at  two  circumferential 
locations.  The  average  clearance  at  design  speed  was  about  0.020 
to  0.025  inches,  approximately  0.6%  of  the  rotor  tip  chord. 

Since  both  sensor  location  and  tip  clearances  were  virtually 
identical  for  both  designs,  these  were  not  considered  as  possible 
explanations  for  any  measured  differences  in  stalling  behavior 
between  the  two  rotors . 

3 . 3  Instriaaentation 

Certain  steady-state  compressor  and  facility  instrumentation 
is  considered  standard  for  all  tests  conducted  in  the  CARL. 

These  were  used  for  the  current  tests  and  will  be  described 
briefly  in  this  section.  More  detail  is  provided  by  Law  and 
Puterbaugh  [29] .  The  main  focus  of  this  section  is  on  the  high- 
response  static  pressure  transducers  used  for  stall  precursor 
detection. 

3.3.1  Test  Facility  Instrumentation 

Rotor  speed  was  measured  by  a  Bentley  Model  306  proximity 
transducer  which  senses  six  grooves  machined  into  the 
gearbox/rotor  driveshaft  coupling.  The  output  was  conditioned 
and  directly  recorded  by  the  data  acquisition  system  (DAS) .  A 
tachometer  provided  visual  indication  of  rotor  speed  accurate  to 
about  10  rpm. 

Inlet  mass  flow  was  metered  through  a  30-inch  venturi  (see 
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Fig.  3.2)  with  a  17.4-inch  throat.  Meter  accuracy  has  been 
calibrated  to  ±0.5%  by  the  manufacturer  (B.I.F.  Industries) . 

Inlet  total  pressure  and  ten^erature  were  measured  just 
downstream  of  the  last  screen  in  the  plenum  (see  Fig.  3.2) . 
Because  of  the  low  plenum  velocities,  compressor  inlet  total 
pressure  was  assvuned  equal  to  plenum  static  pressure  measured  by 
four  static  pressure  taps  manifolded  into  two  pressure  sources. 
At  maximum  flow  rate,  the  error  is  no  worse  than  0.003  psi, 
verified  by  calibration.  Total  temperature  was  averaged  from 
nine  bare  junction  thermocouples  located  in  the  same  axial  plane 
as  the  pressure  taps  at  three  different  radii.  Additionally, 
relative  hiimidity  was  measured  in  the  inlet  for  every  steady- 
state  test  point  and  subsequently  treated  in  the  data  processing 
software  for  the  calculation  of  mass  flow  rate. 

3.3.2  Compressor  InstrusMntation 

Compressor  steady-state  instr\imentation  is  shown 
schematically  in  Fig.  3.9.  A  total  of  276  sensors  were  used  to 
measure  aerodynamic  parameters  at  various  locations  throughout 
the  stage.  Nine  thermocouples  were  mounted  in  the  stator  vane 
leading  edges  auid  eighty  were  located  in  ten  discharge  plane 
rakes .  The  thermocouples  are  of  the  slot  vented  type  shown  in 
Fig.  3.10. 

A  total  of  125  pressures  were  measured  in  the  flowpath; 

36  static  pressures  and  89  total  pressures.  Casing  static 
pressures  were  made  at  fifteen  axial  locations,  twelve  of  these 
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over  the  rotor  tip  (see  Fig.  3.9).  The  total  pressure  probes 
were  all  Kiel  stagnation  tube  design  (Fig.  3.10)  and  were  divided 
in  a  meuiner  identical  to  the  thermocouples.  Both  discharge  plane 
total  pressure  and  temperature  rakes  were  uniformly  distributed 
around  the  circumference  and  spaced  to  divide  a  single  exit  vane 
passage  into  ten  equal  parts . 

The  only  compressor  measurements  used  for  the  current 
analysis  were  the  total  pressures  and  temperatures  measured  by 
the  discharge  rakes  (Fig.  3.9).  These  were  averaged  2uid  used 
with  the  inlet  measurements  in  the  computation  of  stage  pressure 
ratio  and  adiabatic  efficiency  for  each  steady-state  data  point 
acquired  during  the  test  program  (i.e..  Figs.  4.1  and  4.2). 

3.3.3  Calibration  and  Measuranant  Uneartainty 

Since  the  main  concern  of  the  present  work  is  the  analysis 
of  unsteady  pressure  signals^  the  steady-state  pressure  and 
temperature  calibrations  will  not  be  discussed  here.  Details  of 
the  calibrations  are  provided  by  Law  and  Puterbaugh  [29] .  The 
steady-state  measurements  were  used  in  the  current  research  to 
define  where  precisely  on  the  fan  performance  map  the  first  signs 
of  pre-stall  behavior  occur.  Table  3.3  indicates  the  accuracy  to 
which  this  could  be  done.  The  uncertainties  reported  in  Table 

3.3  were  determined  through  a  combination  of  measurement  device 
performeuice  specifications  (for  bias  error)  and  facility 
historical  performance  (for  precision  error) .  Each  measured 
value  (pressure,  ten^erature,  etc.)  was  perturbed  slightly  to 
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obtain  the  performance  sensitivitiea  to  changes  in  various 
parameters . 

3.3.4  Data  Acquisition  System  (DAS) 

Fan  rig  performance  and  calibration  data  were  collected  by 
the  DAS.  The  DAS  consists  of  a  MODCOMP  MODACS  digital  and  analog 
I/O  subsystem,  a  MODCOMP  ATC  commtinications  I/O  subsystem,  and  a 
high  frequency  analog  data  recording  subsystem,  all  controlled  by 
a  host  computer.  The  host  computer  is  a  MODCOMP  Classic  11/15 
16-bit  microprocessor  with  512  kilobytes  of  memory  and  additional 
peripherals.  Additional  information  can  be  foxmd  in  the  report 
by  Law  and  Puterbaugh  [29] . 

Facility  and  rig  controls  and  conqpressor  operating 
parameters  are  all  patched  to  a  control  room  for  monitoring 
during  testing.  A  photograph  of  the  CARL  control  room  is 
provided  in  Fig.  3.11.  The  DAS  software  capabilities  include 
real-time  update  of  con^ressor  performance  parameters,  automated 
data  recording,  and  DAS  health  monitoring. 

3.3.5  Unsteady  Pressure  Measurement 

Eight  high-response  static  pressure  transducers  (Kulite 
model  XCEW-1-187-100  G)  were  located  in  the  conpressor  inlet  for 
detection  of  unsteady  pressure  perturbations.  Table  3.4 
summarizes  the  performance  characteristics  of  the  transducers. 

The  transducers  were  equally  spaced  around  the  circumference 
casing  approximately  0.3  rotor  radii  (2.5  inches)  upstream  of  the 
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rotor  tips  (exact  locations  for  each  rotor  shown  in  Figs.  3.6  and 
3.7)  .  Actual  circiomferential  location  of  each  transducer  is 
indicated  in  Fig.  3.12.  As  noted  in  the  Appendix^  the  equal 
spacing  around  the  annulus  simplifies  the  Fourier  analysis  that 
is  useful  in  determining  the  presence  of  modal  waves . 

Since  the  present  work  was  concerned  with  unsteady  pressure 
signals,  the  Kulites  were  operated  in  a  differential  mode.  The 
reference  side  of  each  transducer  was  manifolded  to  a  static 
pressure  tap  located  at  approximately  the  same  measurement  plane 
as  the  Kulites.  A  long  t\ibe  length  was  used  to  effectively 
filter  any  unsteadiness  from  the  reference  side.  Thus,  measured 
signals  were  true  unsteady  pressure  perturbations. 

The  path  of  each  of  the  Kulite  signals  is  shown  in  Fig. 

3.12.  A  B&F  model  input  conditioner  provided  15-volt  excitation 
and  completed  the  electronic  bridge.  Variable  gain  anqplifiers 
(Bell  and  Howell,  and  Meff)  were  adjusted  to  allow  for  maximum 
sensitivity  without  driving  the  background  noise  levels  too  high. 
The  signals  of  interest  were  expected  to  be  on  the  order  of  one- 
tenth  of  a  psi,  thus  it  was  desirable  to  keep  the  background 
noise  below  this  level.  After  amplification,  the  signals  were 
analog  recorded  (FM  on  100  khz  carrier)  on  a  14 -channel  Bell  and 
Howell  model  VR-3700B  recorder  and  patched  to  oscilloscopes  for 
on-line  monitoring.  A  spectrxim  analyzer  was  also  available  for 
on-line  use.  All  data  processing  was  performed  post -test,  and  is 
described  in  the  next  chapter. 

In-place  calibration  of  the  Kulites  was  accomplished  using  a 
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Druck  DPI  510  pressure  controller/calibrator,  sequentially 
applying  the  calibration  pressure  to  each  transducer.  An^lifier 
gains  were  adjusted  (resulting  gains  were  about  200)  so  that  5.0 
psig  provided  a  2.0-volt  DC  transducer  output.  Static  (non¬ 
rotating)  noise  levels  with  these  an^lification  settings  were  no 
greater  than  10  mV  peak-to-peak  for  each  of  the  eight  channels. 

Each  transducer  was  calibrated  with  5.0,  2.0,  1.0,  0.5,  0.2, 
and  0 . 1  psig  control  pressures .  The  calibration  data  is 
stammarized  in  Table  3.5.  An  inqaortant  consideration  from  the 
calibration  process  was  that  any  signal  less  than  approximately 
0.1  psi  would  be  lost  in  the  noise.  All  the  calibration  curves 
were  linear  down  to  about  0.5  psig.  In  fact,  only  Kulites  2,  3, 
and  5  displayed  significant  non-linearities  (see  Table  3.5)  at 
the  very  low  pressures  (<0.5  psig).  Despite  the  non-linearities 
of  the  three  transducers,  linear  calibration  curves  (all  with  a 
slope  of  2.5)  were  applied  to  each  Kulite  signal  for  the 
following  reasons : 

(1)  The  unfiltered  measured  pressures  all  contained  high 
anqplitude,  high  frequency  components  (i.e.,  blade  passing 
frequency)  that  placed  the  Kulites  in  their  linear  operating 
range  over  most  of  the  low-frequency  period  of  interest  (the 
fundamental  frequency  of  interest  is  approximately  50%  of  the 
rotor  frequency) . 

(2)  The  prime  concern  is  frec[uency  content  rather  than 
exact  aiiq>litude  definition. 
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4.0  TEST  PROCBDDRK  AND  DATA  PROCB88ZNO  OP  UN8TBADY  PRB88I1NB8 


4 . 1  T«st  Proo«dur« 

All  of  the  Kulite  data  were  obtained  at  both  steady-state 
and  trzuisient  operating  conditions  at  five  different  speeds  (60, 
70,  80,  90,  and  100%  design  speed)  for  each  rotor.  For  both 
rotors,  a  detailed  steady-state  performance  map  was  generated 
prior  to  the  acquisition  of  any  unsteady  pressures.  The  maps, 
shown  in  Figs.  4.1  (\inswept.  Rotor  4)  and  4.2  (backward  swept. 
Rotor  6) ,  were  used  for  determining  throttle  settings  for 
starting  the  transients  to  stall  and  for  the  near-stall  steady 
operating  conditions. 

The  steady-state  data  were  obtained  very  systematically. 
Upon  reaching  a  desired  speed,  a  three-minute  settling  time  was 
used  to  allow  system  transients  auid  thermals  to  steady  out. 

After  this,  a  throttle  position  was  set  and  another  two  minutes 
of  settling  time  was  used  before  acquiring  thirty  seconds  of 
Kulite  data.  The  procedure  was  repeated  for  each  change  in 
throttle  position  and  compressor  speed. 

Only  one  throttle  closure  rate  was  available  for  obtaining 
the  compressor  transient  data  to  stall.  At  each  speed,  the 
throttle  was  initially  set  so  that  the  compressor  was  at  some 
nominal  tinchoked  operating  condition  (about  7.0%  greater  flow 
than  that  at  the  near-stall  steady-state  condition)  .  The  FM 
recording  tape  was  then  turned  on  emd  Kulite  data  recorded  while 
the  throttle  was  continuously  closed.  Once  stall  was  detected 
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(both  audibly  and  by  monitoring  th®  Kulitos) ,  th®  surg®  valv®  waa 
immadiataly  op®n®d  and  th®  stall  waa  claarad.  A  t2^ical  throttl® 
closur®  to  stall  (treuisiant)  took  batwaan  thirtaan  and  fiftaan 
saconds.  Tha  throttle  closure  rata  was  approximately  -0.3 
Ibm/s/Sf  calculated  as  defined  below: 

dy/_  ^acall~^iniciml  ^_Q  ^  Ibnt/ 3 

dt  At  '  s 

During  the  transients,  it  was  observed  that  tha  con^rassor 
stalled  at  a  lower  mass  flow  than  the  steady-state  stall  point 
(the  indicated  throttle  position  was  about  5%  more  closed  during 
transient  stall) .  This  was  not  unexpected  given  the  rather  fast 
closure  rate  (compare  eqn.  4.1  with  the  rates  used  by  Gamier 
[28]  discussed  in  section  2.4)  and  the  well-known  fact  that  most 
separation  phenomena  seem  to  have  a  time  constant  associated  with 
them.  Figure  4.3  schematically  demonstrates  the  transient 
behavior  of  both  an  isolated  airfoil  and  conpressor.  This 
unsteady  flow  behavior  helps  explain  why  rotating  disturbances 
could  be  detected  prior  to  stall  during  the  transients,  but  were 
not  discernible  in  the  majority  of  steady-state  points  (discussed 
in  Chapter  5.0) .  Unfortunately,  other  than  the  inlet  Kulites,  no 
transient  measurements  were  availsUble  for  the  test  programs. 
Consequently,  accurate  location  of  the  transient  stall  line  could 
not  be  determined. 

4 . 2  Data  Sasqpling 

As  noted  previously,  all  unsteady  pressure  signals  were 
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recorded  on  analog  tape  and  poat-proceased  after  the  teating  waa 
complete.  The  tapea  were  played  back  on  a  Honeywell  model 
Ninety-aix  Magnetic  Tape  Syatem  and  all  aignala  were  low-paaa 
filtered  uaing  a  Rockland  Multi-channel  Filter  Syatem  816.  The 
filtering  waa  uaed  to  prevent  maaaive  file  atorage  problema,  and 
becauae  the  maximum  frequenciea  of  intereat  were  expected  to  be 
approximately  600-700  Hz  (the  third  harmonic  of  50%  of  the  rotor 
frequency,  the  fundamental  frequency  of  intereat) .  Conaequently, 
an  8 -pole  Butterworth  low-paaa  filter  with  a  cutoff  frequency, 
fc,  of  1  kHz  waa  configured.  Filter  characteriatica  are  ahown  in 
Fig.  4.4.  Alao,  since  the  concern  waa  on  the  preaaure 
perturbationa,  any  DC-biaa  waa  removed  from  the  aignala. 

Taking  the  sanqpling  theorem  (accurate  reproduction  of  an 
analog  signal  requires  swpling  greater  than  twice  the  highest 
frequency  of  intereat,  Stearns  and  Hush  [33])  and  filter 
characteristics  (Fig.  4.4)  into  account,  the  signals  were  san^led 
at  a  rate  of  4  kHz .  Note  that  a  4  kHz  sanqple  rate  places  the 
Nyquist  frequency  (one-half  the  sample  rate)  of  2  kHz  at  twice 
the  filter  cutoff  frequency,  2fc.  Any  signal  greater  than  2  kHz 
and  sampled  at  4  kHz  has  the  potential  to  be  artifically 
introduced  as  a  lower  frequency  (aliasing) ;  however,  according  to 
Fig.  4.4,  a  frequency  at  2fc  is  down  in  amplitude  by  48  db. 

Thus,  the  4  kHz  sample  rate  was  chosen  as  a  good  tradeoff  between 
controlling  file  size  emd  minimizing  any  aliasing.  Examination 
of  the  spectrxims  of  the  unfiltered  and  filtered  signals  revealed 
that  no  aliasing  occured.  Thus,  with  8  chaxmels  digitized  at  4 
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kHz,  sampling  a  10-second  record  resulted  in  a  2.56-Mbyte  data 
file  (8  bytes  per  saunple  per  channel)  .  All  data  were  stored  in 
standard  ASCII  format  with  each  column  containing  the  sampled 
filtered  signal  from  one  of  the  eight  Kulites. 

4 . 3  Additional  Digital  Filtering 

In  addition  to  the  analog  filtering,  unwanted  frequencies 
were  removed  by  using  digital  finite  impulse  response  (FIR) 
filters .  These  were  chosen  for  their  desirable  characteristic  of 
linear  phase  shift  (constant  time  delay) . 

Upon  examination  of  the  raw  analog  filtered  traces  and 
spectrums  of  those  traces,  it  was  found  that  several  unwanted 
frequencies  were  present  from  many  different  sources.  All  traces 
showed  a  1/rev  rotating  wave  (i.e.,  at  the  rotor  frequency)  and  a 
very  strong  low  frequency  planar  wave  (and  several  harmonics) . 
Further,  several  channels  had  strong  60  Hz  (and  harmonics) 
electrical  noise  components .  Figure  4 . 5  shows  an  example  of  a 
frequency  spectrxim  that  was  obtained  from  Kulite  #2  at  90%  speed 
(310  Hz) .  All  of  the  just-mentioned  frequencies  are  evident  in 
the  spectrum.  The  planar  wave  with  fundamental  frequency  at 
approximately  13.5  Hz  is  believed  to  be  a  system  resonauice.  It 
had  a  significeuit  effect  on  the  results  and  is  discussed  further 
in  the  next  chapter. 

Among  the  many  alternatives  for  handling  the  \inwanted 
frequencies,  it  was  decided  to  digitally  band-pass  filter  the 
data  aroiond  the  rotating  stall  frequency.  This  limited  the 
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examination  of  the  modal  waves  to  the  first  harmonic  (SFC  #1) , 
but  was  done  for  the  following  reasons: 

(1)  It  was  the  simplest  means  for  removing  the  majority  of 
the  unwanted  frequencies . 

(2)  Examination  of  the  first  three  rotating  modes  (via 
Fourier  techniques  discussed  in  section  4.4.2)  revealed  that  most 
of  the  modal  wave  energy  appeared  to  be  contained  in  the  first 
harmonic  and  the  additional  filtering  enhanced  its  features. 

(3)  The  frequency  of  the  first  harmonic  for  all  speeds  was 
approximately  50%  of  rotor  frequency  (the  rotating  stall 
frequency) .  Consequently,  the  second  modal  wave  harmonic  showed 
up  right  around  the  1/rev  frequency  making  it  difficult  to 
resolve. 

Shown  in  Fig.  4.6  is  the  same  trace  as  in  Fig.  4.5,  but 
band-pass  filtered  between  130  Hz  and  230  Hz.  Comparison  of  the 
two  figures  reveals  no  attenuation  of  the  frequencies  in  the 
range  of  interest.  Additionally,  no  new  information  has  been 
introduced  as  a  result  of  the  filtering.  The  before  and  after 
spectra  of  each  filter  application  were  checked  in  this  manner  to 
ensure  no  erroneous  results  were  introduced.  Table  4.1  contains 
a  summary  of  the  frequency  band-pass  applied  to  the  data  at  each 
compressor  speed.  In  each  case,  a  248-order,  Hamming-windowed, 
FIR  filter  was  applied  using  MATLAB  [30]  tools.  Unless  otherwise 
indicated,  all  plots  presented  in  the  next  chapter  were  generated 
from  the  band-pass  filtered  data  with  frequency  range  given  in 
Table  4.1. 
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4.4  Forms  of  Data  Prasantation 


Three  basic  forms  of  data  presentation  were  used  to  aid  in 
the  analysis  of  the  measured  unsteady  pressures.  These  included 
plots  of  the  filtered  raw  traces,  magnitude  and  phase  of  spatial 
Fourier  coefficients  (SFCs) ,  and  power  spectral  densities  (PSDs) 
of  the  SFCs .  Figure  4 . 7  summarizes  the  signal  processing 
techniques  used  for  characterizing  the  stalling  process.  MhTLAB 
[30]  was  the  software  package  used  for  the  signal  processing  and 
plot  generation.  Each  form  provided  valuable  information  as  is 
seen  in  the  next  chapter .  The  details  of  each  are  provided 
below . 


4.4.1  Filtered  Raw  Pressure  Signals 

Prior  to  any  Fourier  analysis,  the  filtered  raw  pressure 
signals  were  examined  extensively.  In  particular,  the  last 
20-100  rotor  revolutions  (rotor  revs)  prior  to  stall  were 
carefully  studied  for  any  evidence  of  a  stall  precursor. 
Calibration  curves  were  applied  to  all  pressure  signals  so  that 
plots  of  the  filtered  data  traces  reveal  actual  pressure 
perturbations  in  psi  (rather  than  volts).  Figure  4.8  is  an 
ex2unple  of  a  raw  signal  plot  which  clearly  indicates  the  presence 
of  a  rotating  disturbance.  The  traces  were  artificially 
separated  to  help  view  the  individual  signals  (otherwise  they 
would  lay  on  top  of  each  other)  and  track  the  rotating  phenomena. 
In  all  the  plots  which  display  the  eight  individual  signals, 
Kulite  #1  (see  Fig.  3.12)  provides  the  bottom  trace  and  Kulite  #8 
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the  top  trace. 


4.4.2  Spatial  Fourier  Coefficient* 

The  pressure  signals  were  spatially  resolved  through  the  use 
of  the  discrete  Fourier  transform  (DFT) .  A  review  of  Fourier 
techniques  with  application  to  the  current  work  is  provided  in 
the  Appendix.  By  taking  the  period  to  be  eighty  the  SFCs  can  be 
computed  at  each  discrete  time  from: 

7  -J2*nk 

t)e  *  (4.2) 

°  n«0 

where  k  is  the  mode  number,  or  wave  harmonic.  As  discussed  in 
the  Appendix,  the  c,,'s  are  con^lex  conjugates,  but  for  the  case 
of  eight  equally  spaced  measurements,  only  the  c.,^  of  the  first  3 
modes  (c.^/  c.,,  c.,)  need  be  computed.  Thus,  it  is  seen  from  Eqn. 
(4.2)  that  eight  pressure  traces  are  reduced  to  a  single  coo^lex 
number  at  each  insteuit  in  time.  The  complex  number  (SFC) 
contains  the  desired  phase  and  eunplitude  information  for  each 
spatial  wave  number,  k,  which  ceui  then  be  tracked  in  time. 

If: 

c,9x*jy  (4.3) 

Then: 


I  Cjf  I  +y  *  (4.4a) 

Zc*»ran-M-^)  (4.4b) 

X 
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If  a  sinusoid  is  present,  the  phase,  Lc^,  represents  the  angular 
position  of  its  peak. 

Plotting  the  amplitude  and  phase  of  each  SFC  as  functions  of 
time  helps  reveal  any  rotating  perturbations  that  might  be 
present  in  the  system.  Of  these,  the  unwrapped  phase  plots^  can 
be  the  most  revealing  as  a  rotating  disturbance  propagating  at  a 
constant  speed  will  show  up  as  a  constant  slope.  However,  care 
must  be  taken  in  interpreting  these  plots  as  the  Lc^  will  "lock 
on"  to  the  most  dominant  rotating  disturbance. 

Figure  4 . 9  illustrates  the  usefulness  and  potential  pitfalls 
of  applying  the  spatial  DFT.  The  plots  show  70  rotor  revs  up  to 
and  including  the  start  of  a  rotating  stall  at  90%  speed  (stall  6 
t»464) .  Both  phase  plots  indicate  the  presence  of  a  rotating 
disturbance.  However,  calculations  of  the  slopes  reveal  that  in 
Fig.  4.9a,  the  disturbance  is  propagating  at  the  rotor  frequency 
(l/rev>«310  Hz)  .  This  is  further  evidenced  by  the  change  in  slope 
at  approximately  t«461  when  the  rotating  stall  frequency 
dominates.  In  Fig  4.9b,  the  undesired  rotor  frequency  has  been 
removed  prior  to  performing  the  DFT,  and  the  phase  plot  now 
reveals  the  desired  rotating  stall  frequency  prior  to  stall  (no 
change  in  slope  at  t»461) .  It  is  also  of  interest  to  con^are  the 
aunplitude  plots  of  both  cases.  Note  that  the  1/rev  disturbance 
is  the  major  contributor  to  the  amplitude  prior  to  the  actual 
stall  as  evidenced  by  the  reduced  magnitude  plot  (Fig.  4.9b) 

*Phase  plots  are  typically  confined  between  it  and  -n.  For 
the  case  of  a  rotating  disturbance,  it  is  convenient  to  unwrap 
the  phase  and  allow  it  to  grow  with  time. 
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after  the  1/rev  component  is  removed. 


4.4.3  Power  Spectral  Densities 

By  performing  a  temporal  Fourier  transform  or  standard  fast 
Fourier  transform  (FFT)  on  each  SFC,  the  power  spectral  density 
(PSD)  of  an  SFC  can  be  computed.  Examination  of  the  PSDs  can 
reveal  pre-stall  waves  that  otherwise  would  be  masked  by  other 
rotating  phenomena,  for  example,  a  strong  1/rev  signal.  Figure 
4.10  is  an  example  of  a  PSD  taken  on  the  same  data  used  to 
generate  Fig.  4.9a  (but  not  including  the  stall  itself).  Note 
that  in  addition  to  the  expected  strong  peak  at  1/rev,  a  peak  is 
revealed  at  approximately  50%  rotor  frequency.  No  evidence  of 
this  frequency  was  apparent  in  the  phase  plot  of  Fig.  4.9a. 

When  computing  the  spectrum,  only  data  up  to  the  stall  event 
must  be  used.  If  any  of  the  in-stall  data  is  included,  the  stall 
will  dominate  the  spectrum  and  no  conclusions  will  be  able  to  be 
made  regarding  the  small  amplitude  pre-stall  behavior.  To 
alleviate  this  problem,  PSDs  were  performed  on  the  near-stall 
steady-state  data  (as  close  to  stall  as  possible  without 
inadvertantly  dropping  into  stall)  and  the  portion  of  the 
transient  data  up  to  (but  not  including)  stall.  Also,  as 
described  in  the  Appendix,  because  the  SFCs  are  coxtplex  numbers, 
the  PSDs  are  not  necessarily  symmetric  about  zero  frequency. 

Only  those  PSD  peaks  which  do  not  exist  at  identical  positive  and 
negative  frequencies  indicate  rotating  disturbauices .  Thus,  both 
the  positive  euid  negative  frequencies  were  overplotted. 
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5.0  RESULTS  AND  DISCUSSION 


5.1  Introductory  Rwaarks 

Presentation  of  results  is  provided  primarily  in  sections 
5.3  and  5.4  along  with  some  discussion.  Because  of  their  effect 
on  the  results r  the  low  frequency  planar  oscillations  mentioned 
in  section  4.3  are  discussed  first  in  the  next  section. 

Additional  discussions  regarding  the  effects  of  compressor 
speedline  characteristic  and  shock  waves  are  presented  in 
sections  5.5  and  5.6.  In  the  last  section,  results  are  compared 
with  those  from  previous  investigations. 

5.2  Low  Frequency  Planer  Oscillations 

As  noted  in  section  4.3,  planar  waves  with  fundamental 
frequency  of  about  13.5  Hz  were  evident  in  all  transient  and 
near-stall  steady  data  obtained  on  both  rotors.  As  exemplified 
in  Fig.  4.5,  the  second  harmonic  of  this  wave  (27  Hz)  was 
dominant  at  all  rotor  speeds  except  100%  speed  for  both  designs. 
At  100%  speed,  the  fundamental  frequency  dominated.  Although  the 
precise  nature  of  the  planar  oscillations  could  not  be  determined 
(these  type  oscillations  have  been  observed  in  previous 
investigations;  see  Greitzer  [31]),  they  are  believed  to  be  the 
result  of  some  overall  system  resonance  (rig  and  facility)  for 
the  following  reasons: 

(1)  the  frequency  did  not  change  with  compressor  speed. 

(2)  the  waves  were  seen  in  the  near-stall  steady-state 
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data^  thus  eliminating  the  notion  that  they  might  be  induced  by 
the  throttle  ramp  rate. 

(3)  the  facility  is  closed-loop  and  believed  to  be  more 
susceptible  to  resonances  than  an  open-loop  facility. 

In  order  to  verify  the  system  resonance  theory,  the 
compressor  zmd  inlet  and  exit  ducting  were  modeled  as  a  Helmholtz 
resonator  in  a  manner  similar  to  Emmons,  et  al.  [21]  emd  Greitzer 
[31]  .  The  Helmholtz  frequency,  o),  is  defined  as: 


Ci)®a. 


N  V-c 


(5.1) 


where  and  L«,  are  the  area  euid  length  of  the  equivalent 
compressor  with  constant  area  inlet  duct,  is  the  discharge 
plenxim  volume,  and  a  is  the  speed  of  sound.  For  the  current  set 
of  tests,  the  discharge  volume  was  fixed  and  the  throttle  choked; 
thus,  Vp  was  calculated  with  a  high  degree  of  confidence. 

However,  choosing  euid  were  somewhat  arbitrary,  especially 
with  a  closed  loop  facility.  Two  attempts  were  made,  first 
including  inlet  volume  only  up  to  the  screens  in  the  inlet  plenvim 
(see  Fig.  3.2),  then  adding  the  remaining  inlet  plenum  volume  up 
to  the  perforated  cone  (Fig.  3.2) .  The  first  attempt  produced  an 
(0  *149  Hz,  the  second,  (D  ">80  Hz.  The  question  as  to  which  inlet 
volume  should  or  should  not  be  included  is  undecided,  and 
certainly  by  including  enough,  a  13.5  Hz  (or  27  Hz)  resonator 
could  be  established. 

The  presence  of  the  planar  oscillations  con^licated  the 
signal  processing  procedures  used  for  detection  of  the  rotating 
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phenomena  of  interest  here.  This  complication  arises  because  the 
oscillations  were  observed  close  to  the  onset  of  stall  and  were 
consistently  detected  before  any  rotating  perturbations  were 
observed  (i.e.^  platnar  waves  could  be  seen  in  the  near-stall 
steady  data,  rotating  waves  were  not) . 

As  seen  in  Fig.  5.1a,  the  planar  waves  were  initially  of 
small  amplitude  and  grew  as  the  machine  was  continually  throttled 
to  stall.  Even  when  the  data  were  high-pass  filtered  above 
60  Hz,  the  effect  of  the  waves  was  clear  (Fig.  5.1b);  the 
oscillations  effectively  throttled  the  stage  into  and  out  of 
stall.  A  clearer  picture  of  this  is  provided  by  Fig.  5.2,  which 
shows  the  in-stall  behavior  at  80%  speed.  The  behavior  shown  is 
quite  consistent  with  surge-like  oscillations  (see,  for  exai^ple, 
Greitzer  [31]).  Of  particular  interest  to  the  current  work  is 
the  small -amplitude  "surging"  that  occurred  prior  to  stall  and 
its  influence  on  the  rotating  disturbances.  This  influence  is 
shown  schematically  in  Fig.  5.3  and  is  evidenced  in  the  majority 
of  the  results  presented  in  the  next  section. 

5.3  Analysis  of  Results 

Because  the  rotors  studied  are  so  close  in  design,  it  is 
appropriate  to  present  results  via  speedline  (for  both  rotors) 
rather  than  by  rotor.  This  allows  for  direct  comparison  between 
designs  and  facilitates  a  discussion  regarding  any  significant 
differences  in  stalling  behavior.  Signal  processing  of  the  data 
was  performed  as  described  in  Chapter  4.0.  Some  exan^les  of 
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signal  traces  without  the  additional  filtering  (Table  4.1)  are 
also  provided.  Plots  presented  as  fxonctions  of  time  use  rotor 
revolutions  (rotor  revs)  as  the  time  scale  since  this  is  the 
scale  of  the  phenomena  of  interest.  In  most  time  plots ^  only  the 
last  1-2  seconds  prior  to  the  stall  event  are  shown  since  the 
current  investigation  is  concerned  with  stall  inception. 

Rotating  stall  frequencies  for  all  cases  were  between  50-56%  of 
rotor  frequency.  A  single  stall  cell  was  identified  in  all  the 
cases  exaunined.  As  noted  in  section  1.2,  both  progressive  and 
abrupt  rotating  stall  were  observed.  Progressive  stall  was  seen 
at  low-speed  operation  (60%  speed  for  the  unswept  rotor,  #4,  and 
60  and  70%  speed  for  the  swept  rotor,  #6) ,  as  would  be  expected 
(Greitzer  [12]). 

At  all  speeds  for  both  rotors,  stall  was  preceded  by  small 
amplitude  rotating  waves  which  grew  rapidly  once  the  onset  of  the 
instaQ^ility  was  reached.  Thus,  the  plots  presented  for  either 
rotor  are  meant  to  be  representative  of  the  stall  inception 
behavior  of  both  designs.  The  only  significant  difference  was  in 
the  eunount  of  time  the  modal  waves  could  be  detected  prior  to 
stall .  The  differences  in  stall  warning  time  will  be  discussed 
qualitatively  in  this  section,  and  more  quantitatively  in  the 
following  section  (section  5.4) . 

The  data  digitizing  system  described  in  section  4.2  became 
inoperative  before  any  steady-state  data  for  Rotor  6  could  be 
sampled.  In  this  sense,  the  data  emalysis  was  not  as  conq>lete  as 
it  could  have  been.  However,  given  the  data  that  were  examined 


46 


(Rotor  4  steady-state  and  tremsient  and  Rotor  6  transient) ,  a 
high  level  of  confidence  is  placed  in  the  results  and  conclusions 
that  follow. 

5.3.1  60%  Design  Speed 

This  was  the  only  speed  where  the  relative  rotor  flow 
velocities  were  subsonic  across  the  entire  span.  As  noted  above, 
a  progressive  rotating  stall  was  observed  for  both  designs  at 
this  speed.  The  data  traces  for  Rotor  4  (\jnswept)  were 
previously  presented  in  Fig.  5.1;  the  traces  for  Rotor  6  were 
quite  similar  and  are  not  shown.  The  progressive  nature  of  the 
stall  cam  be  seen  in  Fig.  5.1b  as  the  pressure  perturbations 
experience  a  relatively  slow  growth  (compare  the  traces  in  Fig. 
5.1b  to  those  of  Fig.  5.7b).  As  indicated  in  Fig.  5.4,  initial 
signs  of  the  stall  became  evident  at  around  t*196.  Note  the 
small  an^litude  rotating  wave  (<•0.2  psi  peak-to-peak)  at 
approximately  50%  rotor  frequency  (0.5/rev)  emerging  at  t»196  and 
being  attenuated  by  t«201. 

From  a  stall  warning  viewpoint,  this  case  of  a  progressive 
stall  is  relatively  uninteresting  as  one  cem  throttle  in  and  out 
of  stall  without  sudden  cheuiges  in  compressor  performance. 
However,  for  completeness,  the  amplitude  and  phase  plots  for  SFC 
#1  are  provided  in  Fig.  5.5  and  do  indicate  the  presence  of  the 
rotating  disturbance  from  about  t»190  on.  Note  that  traces  of 
the  disturbance  can  be  seen  prior  to  t*190. 

Figure  5.6  shows  the  PSD  of  SFC  #1  at  the  near-stall  steady 
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condition  (DV053)  and  next  higher  flow  condition  (DV050)  for 
which  data  were  obtained.  The  insert  identifies  the  compressor 
operating  conditions  (%  design  pressure  ratio  versus  %  design 
flow)  corresponding  to  these  two  throttle  settings.  The  peak  at 
110  Hz  (about  51%  of  rotor  frequency;  see  Fig.  5.5  phase  plot)  at 
the  near-stall  condition  provides  further  evidence  of  the  pre¬ 
stall  rotating  waves.  No  discernible  peak  at  110  Hz  is  evident 
at  throttle  setting  DV050;  hence,  no  rotating  wave  behavior  is 
apparent.  This  was  the  only  speed  in  which  a  pre-stall 
indication  was  evident  in  the  near-stall  steady-state  data. 

5.3.2  70%  Design  Speed 

As  evidenced  in  Fig.  5.7,  the  stall  at  70%  speed  for  Rotor  4 
was  abrupt  and  relatively  clean  (the  effects  of  the  low  frequency 
oscillations  were  minimal) .  The  planar  oscillations  described  in 
section  5.2  were  still  present  (their  effect  can  be  seen  best  in 
Fig.  5.7b  through  low  frequency  modulation  of  the  rotating  stall 
frequency) ,  but  were  not  as  strong.  Examination  of  the  last 
twenty  rotor  revs  prior  to  stall  (Fig.  5.8)  reveals  the  stall 
cell  emerging  smoothly  from  small  amplitude  rotating  waves  (at 
about  t*483) .  Prior  to  t-483,  it  is  difficult  to  detect  any 
coherent  rotating  disturbance  in  the  raw  traces. 

Because  the  Kulite  traces  were  relatively  clean  of  \uiwanted 
frequencies,  higher  order  Fourier  modes  were  also  examined  at 
this  speed.  Shown  in  Fig.  5.9  are  the  an^litude  and  phase  plots 
of  the  first  two  spatial  modes  (SFC  #3  contained  virtually  no 
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energy)  .  As  indicated  in  the  magnitude  plot  prior  to  stall 
(stall  at  t*485’)  ,  the  energy  content  of  the  first  two  modes  is 
quite  similar,  with  SFC  #1  displaying  slightly  more  activity. 

Both  phase  plots  appear  to  show  a  reduced  slope  change  at  about 
t»470  (from  1/rev  to  »0. 52/rev) .  The  slope  from  t»470  on 
corresponds  to  a  rotating  disturbance  traveling  around  the 
annulus  at  about  52%  of  the  rotor  frequency  (note  the  slope 
remains  the  same  subsequent  to  stall  onset  at  t*485) ,  indicative 
of  the  impending  stall.  Using  the  phase  plot  data  of  Fig.  5.9, 
it  appears  that  minimal  stall  warning  time  (15  rotor  revs)  is 
provided  at  this  speed  for  Rotor  4. 

Plotting  the  PSD  of  the  spatial  Fourier  coefficients 
provides  another  tool  for  examining  the  stall  inception  behavior 
of  the  compressor.  The  data  obtained  during  the  transient  to 
stall  can  be  used  provided  the  stall  event  itself  is  not 
included.  The  power  spectra  will  then  reveal  any  small  amplitude 
resonances  in  the  system  even  if  they  are  present  only  for  a 
short  duration.  This  was  the  procedure  followed  since  the 
steady-state  data  provided  no  insights  into  the  stalling 
mechanism. 

Shown  in  Fig.  5.10  are  PSDs  of  SFC  #1  at  70%  speed  obtained 
using  the  transient  data  including  stall  (5.10a)  and  up  to  stall 
(5.10b) .  The  same  plots  are  shown  for  the  second  mode  in  Fig. 
5.11.  Because  the  SFCs  are  complex,  their  spectra  are  not 

^or  consistency  at  all  speeds,  the  onset  of  stall  was  taken 
to  be  when  the  amplitude  of  SFC  #1  reached  50%  of  its  final 
value . 
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necessarily  symmetric  (explained  in  the  Appendix  aund  Chapter 
4.0);  consequently,  positive  and  negative  frequencies  are 
overplotted.  Peaks  that  exist  at  identical  positive  and  negative 
frequencies  indicate  a  standing  wave  (note  the  twin  peaks  at  120 
Hz  in  Figs.  5.10b  and  5.11  which  imply  an  amplitude  oscillation 
of  a  standing  wave;  due,  in  this  case,  to  electrical  noise  which 
was  strongly  present  on  some  channels  and  not  evident  at  all  on 
others) .  Thus,  only  those  signals  which  do  not  occur  at 
identical  peaks  indicate  rotating  disturbances.  A  further 
complication  for  this  particular  speed  is  that  the  rotating  stall 
frequency  is  quite  close  to  the  120  Hz  noise  frequency. 
Nevertheless,  the  stall  frequency  did  show  up  to  be  slightly 
higher  (123  Hz)  in  the  data  which  included  the  stall  (Fig. 

5.10a) .  The  same  is  true  in  Fig.  5.11a  for  the  second  harmonic 
(246  Hz) .  Examination  of  Figs.  5.10b  and  5.11b  reveals  no  peaks 
at  the  stall  frequency.  Given  the  120  Hz  complication  and  the 
complexity  added  by  the  low  frequency  plemar  waves,  the  spectra 
offered  no  additional  information  regarding  the  stalling  process . 

The  stall  of  Rotor  6  (swept)  at  70%  speed  was  quite 
different  from  that  of  Rotor  4.  The  behavior  of  Rotor  6  was 
characterized  by  a  progressive  stall  very  similar  to  that 
described  in  section  5.3.1  for  both  rotors  at  60%  speed.  The 
progressive  nature  of  the  stall  is  represented  in  the  SFC  #1 
euiqplitude  plot  of  Fig.  5.12,  where  one  can  see  a  relatively 
gradual  increase  in  the  energy  content.  Also,  the  effects  of  the 
low  frequency  pletnar  oscillations  are  apparent  in  the  anplitude 
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plot  of  Fig.  5.12,  as  can  b«  seen  clearly  from  about  t-570  on 
(note  the  modulations  at  approximately  27  Hz)  .  Conqparison  with 
the  SFC  #1  aunplitude  plot  in  Fig.  5.9  shows  that  Rotor  4  stalls 
quite  eJsruptly  relative  to  the  Rotor  6  stall  at  this  speed,  and 
the  effect  of  the  planar  oscillations  is  minimal.  Further 
evidence  of  a  slowly  developing  (progressive)  stall  is  provided 
in  the  phase  plot  of  Fig.  5.12,  which  shows  a  relatively  constant 
slope  (for  the  entire  200  rotor  revs  plotted)  corresponding  to  a 
rotating  disturbance  at  ed^out  0.5/rev.  Despite  the  differences 
in  the  type  of  rotating  stall  encountered  (abrupt.  Rotor  4; 
progressive.  Rotor  6)  by  the  two  designs  at  70%  speed,  the  basic 
mechaoiism  of  the  stall  development  process  of  the  two  rotors  was 
identical  -  growth  of  small  eunplitude  rotating  waves  (Figs.  5.9 
and  5.12) . 

5.3.3  80,  90,  and  100%  Design  Speed 

At  the  higher  compressor  speeds,  abrupt  stall  (sudden, 
discontinuous  decrease  in  performance  at  stall  onset)  was 
observed  for  both  designs .  As  will  be  shown  in  this  section,  the 
stall  inception  process  was  governed  by  the  growth  of  small 
euaplitude  rotating  waves .  Again,  the  plots  presented  for  either 
rotor  are  meant  to  be  representative  of  the  stall  inception 
behavior  of  both  designs .  It  is  worth  noting  that  evidence  of 
the  upcoming  stall  could  be  seen  in  the  unfiltered  (no  analog  or 
digital  filtering)  raw  Kulite  signals.  At  the  higher  rotational 
speeds,  the  pressure  traces  were  dominated  by  a  strong  signal  at 
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the  blade  passing  frequency  (believed  to  be  caused  by  shock  waves 
passing  the  Kulites) .  As  stall  was  approached,  the  1/blade 
signal  started  oscillating  in  both  frequency  and  amplitude, 
approximately  1-2  seconds  prior  to  stall  onset.  This  was 
observed  both  on-line  while  testing  and  during  playback  of  the 
data.  This  observation  is  discussed  further  in  section  5.7. 

The  filtered  data  traces  obtained  from  the  80%  speed 
transient  for  Rotor  4  (unswept)  are  shown  in  Figs.  5.13  and  5.14. 
The  longer  period  indicated  between  the  first  two  "pulses”  of 
rotating  stall  (seen  best  in  Fig.  5.2)  was  repeatable  with  the 
S2une  rotor  and  was  also  seen  in  the  Rotor  6  data  at  the  same 
speed.  Development  of  the  rotating  stall  from  the  small 
amplitude  modal  waves  can  be  seen  clearly  in  Fig.  5.14.  Figure 
5 . 15  presents  the  amplitude  and  phase  data  of  the  first  mode  and 
points  out  some  interesting  features.  Using  the  information 
presented  in  Fig.  5.15,  there  does  not  appear  to  be  a  lot  of 
warning  time  prior  to  the  instability  as  evidenced  in  both  the 
magnitude  and  phase  plots.  Note  the  apparent  exponential  growth 
of  the  rotating  waves  (from  t»464  on)  just  prior  to  stall  onset. 
This  is  predicted  by  the  stability  theory  discussed  in  section 
2.3.1  (Gamier  [28]).  The  1/rev  disturbance  indicated  between 
t-423  and  435  is  believed  to  be  eui  anomaly,  as  the  data  was 
filtered  using  a  248-pole  FIR  filter  between  125  and  230  Hz. 

The  spectra  obtained  using  the  transient  data  set  (up  to, 
but  not  including  stall)  for  Rotor  4  at  80%  speed  is  shown  in 
Fig.  5.16.  As  evidenced  by  the  lack  of  any  definitive  peak,  the 
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PSD  does  not  appear  to  offer  any  additional  insights  into  the 
process  leading  up  to  stall.  Several  of  the  peaks  (i.e.,  150  Hz, 
175  Hz)  are  believed  to  be  harmonics  of  the  strong  low  frequency 
planar  oscillations  (fundamental  frequency  of  13.5  Hz)  evident  in 
all  the  data  (see  Fig.  4.5). 

The  stall  behavior  at  90%  speed  provided  one  of  the  more 
interesting  comparisons  between  the  two  designs.  Indeed,  the 
development  of  stall  from  rotating  waves  for  both  rotors  was 
again  observed.  Figures  5.17  and  5.18  are  representative  of  the 
process  for  both  designs.  Small  amplitude  wave  growth  cam  be 
seen  clearly  in  Fig.  5.18.  The  most  interesting  feature  of  the 
stalling  process  at  this  speed  for  both  designs  was  that  the 
waves  were  detected  much  sooner  than  those  for  any  other  speed 
examined.  Evidence  of  this  is  provided  in  Figs.  5.19  (Rotor  4, 
unswept)  and  5.20  (Rotor  6,  swept) .  Both  phase  plots  show 
relatively  long  periods  of  constant  slope  prior  to  stall, 
corresponding  to  a  rotating  disturbance  at  aU?out  50%  rotor 
frequency.  Of  particular  interest  are  the  SFC  #1  magnitude  and 
phase  plots  for  Rotor  6  (Fig.  5.20) .  Both  plots  in  Fig.  5.20 
reveal  significant  energy  contained  in  the  first  spatial  mode  for 
at  least  260  rotor  revs  (0.84  s)  prior  to  stall.  The  reason  for 
this  advanced  stall  warning  time  is  strongly  suggested  by 
examination  of  the  compressor  overall  characteristic  at  90%  speed 
(Fig.  4.2) .  Note  that  the  90%  characteristics  for  both  rotors 
(Figs.  4.1  and  4.2)  exhibit  positively- sloped  regions,  with  the 
speedline  for  Rotor  6  being  most  pronounced.  The  effect  of 
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speedline  shape  on  stall  precursor  behavior  is  discussed  in  more 
detail  in  section  5.5. 

The  PSD  plot  of  SFC  #1  for  Rotor  6  (swept)  using  the 
transient  data  up  to  stall  is  shown  in  Fig.  5.21.  As  with  80% 
speed  (Fig.  5.16),  the  spectra  does  not  reveal  any  distinct  peaks 
that  might  be  indicative  of  small  amplitude  rotating  waves.  The 
lack  of  information  from  Fig.  5.21  does  not  invalidate  the 
results  just  described  from  the  phase  plot  of  Fig.  5.20.  Figure 
5.21  simply  shows  that  the  spectra  of  the  first  Fourier  mode  is 
complicated  and  contains  a  lot  of  frequencies  (rotating,  planar) 
of  approximately  the  same  magnitude. 

At  design  speed,  the  stalling  process  was  identical  to  that 
described  at  the  other  high  speeds  **  abrupt  stall  characterized 
by  growth  of  small  amplitude  circximferentially  travelling  waves. 
This  is  evidenced  in  Figs.  5.22-5.24  for  the  swept  rotor  (Rotor 
€) .  The  PSD  plots  of  SFC  #1  for  Rotor  6  are  shown  in  Fig.  5.25 
using  100%  speed  transient  data  up  to  stall,  and  0.25  and  0.5 
seconds  prior  to  stall .  A  very  distinct  rotating  disturbance  can 
be  seen  at  approximately  170  Hz  in  all  three  plots  (Fig.  5.25) . 
Note  how  the  peak  attenuates  and  widens  as  data  further  from  the 
stall  was  used  (Figs .  5 . 25b  and  c) .  This  attenuation  suggests 
that  the  rotating  waves  are  being  danced  out  as  the  compressor 
moves  away  from  the  stall  point,  consistent  with  the  stability 
theory  described  in  section  2.3.1.  Although  not  shown,  the 
spectra  of  SFC  #1  for  Rotor  4  at  design  speed  was  similar  to 
those  shown  in  Fig.  5.16  and  5.21;  no  definitive  peak  was  seen. 
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5.4  Stall  Warning  Tima 

In  order  to  arrive  at  a  more  quantitative  description  of  the 
stall  warning  time,  the  following  procedure  was  adopted.  The 
phase  data  of  SFC  #1  were  smoothed  using  a  cubic  spline  smoothing 
function  (i.e..  Fig.  5.26).  The  choice  of  the  smoothing  fxjnction 
inevitaUDly  affects  the  stall  warning  time,  but  as  noted  below, 
the  criterion  was  developed  for  comparison  purposes  (the  same 
order  smoothing  function  was  applied  to  each  data  set) .  Since 
the  slope  of  the  phase  plot  corresponds  to  the  frequency  of  a 
rotating  disturbance,  the  slope  (derivative)  of  the  smoothed  dats 
was  calculated  and  plotted  as  a  function  of  time  up  to  and 
including  stall.  The  stall  warning  time  was  then  taken  to  be  the 
last  stretch  of  slope  up  to  stall  that  fell  between  30-70%  rotor 
frequency  (2-4.4  rad/ rev ) .  This  criterion  was  developed  merely 
to  provide  a  consistent  means  of  determining  a  stall  warning  time 
for  comparison  purposes.  In  the  plots  that  follow,  the  criterion 
band  and  the  onset  of  fully-developed  stall  are  indicated,  as 
well  as  the  start  time  taken  for  stall  warning  (arrow) . 

Plots  of  the  phase  slope  of  SFC  #1  for  both  rotors  are  shown 
in  Figs.  5.27-5.34  for  70,  80,  90,  and  100%  design  speed.  The 
€0%  speed  data  were  not  included  because  of  the  previously 
described  progressive  stall  exhibited  at  this  speed.  The 
speedline  characteristics  provided  as  inserts  in  each  figure  help 
stimulate  discussion  in  the  next  section.  Bn^hasis  is  placed  on 
the  shape  of  the  speedlines,  hence  they  are  plotted  as  %  design 
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pressure  ratio  versus  %  design  mass  flow  (refer  to  Figs.  4.1  and 
4.2)  . 

Tzdsle  5 . 1  svunmarizes  the  stall  warning  times  as  calculated 
following  the  ahovB  criterion.  As  can  be  seen  from  the  plots  and 
Table  5.1,  Rotor  6  consistently  provided  more  stall  warning  time 
than  did  Rotor  4.  With  the  exception  of  90%  speed  (Figs.  5.31 
and  5.32),  the  shape  of  the  speedlines  for  each  rotor  were  quite 
similar.  The  advanced  stall  warning  times  for  the  swept  rotor 
(Rotor  6)  suggests  a  potential  effect  of  shock  waves  on  stall 
precursor  waves,  since  the  leading  edge  of  Rotor  6  was  swept  back 
deliberately  to  reduce  shock  wave  strength.  Shock  wave  effects 
are  also  discussed  later  in  section  5.6.  In  addition, 
essentially  all  of  the  plots  show  evidence  of  the  impending  stall 
200-300  rotor  revs  prior  to  the  instability  (i.e.,  slope  within 
the  criterion  for  a  30-60  rotor  rev  period  well  before  the  stall, 
but  followed  by  a  significant  stretch  of  the  slope  outside  the 
criterion) .  Very  good  examples  of  this  can  be  seen  in  Figs.  5.33 
and  5.34  at  design  speed  for  both  rotors. 

5.5  Effect  of  Speedline  Characteristic  on  Stall  Precursor 

Figures  5.35  and  5.36  show  the  steady-state  performance 
characteristics  of  the  two  designs  as  determined  using  the  steady 
conditions  for  which  Rulite  data  were  obtained.  The  high-flow 
condition  (wide-open  throttle)  data  points  are  included  to 
provide  a  better  determination  of  the  speedline  shape.  The  low- 
flow  points  at  each  speed  provide  a  good  indication  of  the 
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steady-state  stall  line.  Also  indicated  is  a  best  guess  for  the 
transient  stall  line  (refer  to  section  4.1)  determined  by 
extending  the  general  shape  of  the  speedlines  as  shown  (better 
indications  of  the  speedline  shapes  ceui  be  seen  in  Figs.  4.1  and 
4.2) .  Examination  of  the  extended  speedlines  offers  a  potential 
explanation  as  to  why  the  small  amplitude  waves  were  not  detected 
at  the  near-stall  steady  conditions,  but  were  detected  in  the 
transients  very  close  to  stall.  By  operating  transiently  to  the 
left  of  the  steady  stall  line,  more  of  the  positively  sloped 
region  of  the  speedlines  is  reached.  It  appears  that  a 
positively  sloped  characteristic  (also  discussed  in  section  2.4) 
could  be  an  important  consideration  for  advamced  pre-stall 
indications . 

This  hypothesis  is  further  substantiated  by  examining  the 
stall  warning  times  (Table  5.1  and  Figs.  5.27-5.34)  in  relation 
to  the  speedline  shapes.  As  seen  in  Figs.  5.29  emd  5.30,  the 
slope  of  the  80%  speedline  for  both  rotors  is  nearly  flat  (maybe 
slightly  positive)  and  the  stall  warning  time  is  minimal  compared 
to  the  other  speeds.  In  sharp  contrast  to  the  80%  speedline 
shape  (and  all  other  speeds  for  that  matter)  is  the  very  large 
region  of  positively  sloped  characteristic  at  90%  speed  for  the 
swept  design  (Rotor  6) .  This  condition  provided  more  than  2.0 
seconds  of  stall  warning  time.  Note  that  from  a  design 
standpoint,  operating  in  a  region  of  large  positive  slope  is 
undesirable  -  performance  penalties  are  large  (see  for  exan^le 
the  efficiency  curve  of  Rotor  6  at  90%  speed.  Fig.  4.2) . 
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However,  operating  at  peak  efficiency  in  a  region  of  flat  or 
slightly  positive  slope  is  desirable  (i.e.,  maximum  pressure 
rise) ,  and  may  be  beneficial  from  a  stall  warning  standpoint . 

Figure  5.37  is  an  attempt  to  provide  some  quantitative 
relationship  between  speedline  slope  and  stall  warning  time.  The 
slopes  were  calculated  by  extending  a  straight  line  through  the 
near-stall  steady-state  data  point  following  the  general  shape  of 
the  speedline.  Results  of  this  rather  crude  approximation  (Fig. 
5.37)  point  out  the  difficulty  encountered  in  trying  to  correlate 
stall  precursor  behavior  with  various  compressor  parameters . 
Another  potential  problem  here  is  that  the  statistical  data 
sample  is  quite  limited  (i.e.,  each  rotor  was  stalled  transiently 
only  once  at  each  speed  tested)  .  Since  the  stalling  process  has 
a  lot  of  random  phenomena  affecting  its  behavior,  it  might  be 
more  appropriate  to  look  at  some  average  behavior  for  many  stalls 
repeated  at  each  speed.  Additionally,  speedline  slope  is  just 
one  factor.  Other  important  considerations  were  discussed  in 
section  2 . 4  and  include  sensor  location  and  sensitivity,  throttle 
transient  rate,  flow  field  distortion,  shock  waves,  and  others. 

5.6  Effect  of  Shock  Have*  on  Stall  Precursor 

Figure  5.38  schematically  shows  the  shock  structure  at 
operation  close  to  stall  of  the  transonic  rotors  tested  in  the 
current  experiments.  At  a  high-flow  operating  condition 
(unthrottled) ,  the  leading  edge  shock  is  attached  and  oblique  to 
the  inlet  relative  flow.  As  the  flow  is  throttled  to  stall,  the 
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shock  becomes  more  normal  to  the  relative  flow  and  moves  upstream 
to  form  a  bow  shock  (Fig.  5.38) .  Also  shown  in  Fig.  5.38  is  the 
location  of  the  Kulites  used  for  stall  precursor  detection  in  the 
current  investigation.  With  this  arrzmgement,  the  modal  waves 
rotating  at  about  0.5  rotor  frequency  must  travel  upstream 
(decaying  exponentially  as  predicted  by  theory.  Gamier  [28])  and 
pass  through  the  shock  structure  rotating  at  the  blade  passing 
frequency  (40  times  the  modal  wave  frequency)  in  order  to  be 
detected.  Within  this  complicated  flow  field,  a  quantitative 
description  regarding  the  effect  of  shock  waves  is  difficult; 
however,  qualitatively,  attenuation  of  the  rotating  waves  would 
appear  reasonable. 

As  indicated  in  Table  5.1,  Rotor  6  consistently  displayed 
more  stall  warning  time  than  did  Rotor  4,  even  in  cases  when  the 
speedline  shapes  were  virtually  identical.  Figure  5.39  shows  the 
performance  characteristics  of  both  rotors  on  the  same  map  to 
provide  a  good  comparison  of  the  speedline  shapes.  A  good 
example  can  seen  in  the  80%  speed  data  for  both  designs;  Rotor  6 
had  more  than  double  the  warning  time  of  Rotor  4,  yet  the 
speedline  shapes  were  quite  similar.  In  fact,  all  the 
characteristics  were  very  similar  in  shape  except  at  D0%  speed. 
Also,  the  tip  clearances  and  sensor  location  for  both  designs 
were  essentially  the  same,  as  noted  in  section  3.2.2.  Thus,  it 
appears  that  the  reduced  shock  strength  of  the  swept  rotor  (Rotor 
6)  had  a  positive  effect  on  the  ability  to  detect  the  pre-stall 
waves  and  supports  the  discussion  in  the  preceding  paragraph. 


59 


Although  not  the  focus  of  this  research,  it  is  also 
interesting  to  note  that  the  swept  rotor  stalls  appeared  to  be 
less  severe  thaui  the  unswept  rotor  stalls.  At  70%  speed.  Rotor  6 
exhibited  a  progressive  stall  (Fig.  5.12),  while  the  Rotor  4 
stall  was  more  abrupt  (Fig.  5.9) .  Additionally,  the  SFC  #1 
amplitude  plots  for  the  swept  rotor  showed  consistently  less 
energy  in  fully-developed  stall  than  did  the  \inswept  design  (for 
exaunple,  see  Figs.  5.19  and  5.20).  The  reduced  severity  of  the 
Rotor  6  stalls  appears  to  confirm  the  weakened  shock  structure  of 
the  swept  design. 

5.7  Comparison  With  Previous  Investigations 

Compressor  stall  inception  for  the  two  rotors  examined  was 
characterized  by  growth  of  small  amplitude  rotating  waves.  As 
measured  0.3  rotor  radii  upstream  at  the  blade  tips,  the 
perturbations  were  initially  0.2  psi  peak-to-peak .  The  waves  did 
not  grow  significantly  until  just  prior  to  the  instability,  when 
exponential  growth  into  fully-developed  stall  occurred  very 
rapidly,  typically  within  6-10  rotor  revolutions. 

The  features  described  in  the  preceding  paragraph  support 
the  modal  wave  representation  of  the  stall  inception  process, 
schematically  shown  in  Fig.  2.3.  Results  are  q[ualitatively  quite 
similar  to  those  obtained  from  recent  experiments  on  low-speed 
compressors  (McDougall,  et  al.  [2]  and  Gamier,  et  al.  [3])  and 
one  high-speed  compressor  (Gamier,  et  al.  [3]).  Figure  2.4 
provides  an  exeunple  of  Gamier' s,  et  al.  [3]  results  from  a  low- 
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speed  machine.  The  similarities  of  the  results  from  the  current 
work  (see,  for  example  Figs.  5.19  and  5.20)  with  the  features 
shown  in  Fig.  2.4  are  readily  apparent.  The  planar  oscillations 
evident  in  the  present  results  probably  affected  when  the  modal 
waves  could  initially  be  detected  (see  Fig.  5.3) .  Regardless, 
the  basic  mechanism  of  stall  inception  was  clear  ^  stall  evolved 
smoothly  from  small  eimplitude  circumferentially  travelling  waves. 

It  is  worth  mentioning  that  the  high-speed  results  presented 
by  Gamier,  et  al.  [3]  showed  the  second  harmonic  of  the  rotating 
waves  (SFC  #2)  to  be  at  least  as  important  as  the  first  mode. 

The  present  work  focused  on  the  first  mode  for  reasons  discussed 
in  section  4.3  (a  lot  of  unwanted  frequencies,  most  of  the  energy 
appeared  to  be  in  SFC  #1,  etc.).  The  higher  harmonics  were 
exeunined  only  briefly  and  probably  warrant  more  study.  Given 
that  the  rotating  stall  frequencies  of  both  designs  were 
approximately  50%  of  rotor  frequency,  exeuaination  of  higher  order 
modes  would  require  additional  signal  processing  in  order  to 
resolve  the  second  harmonic  and  the  1/rev  signal.  Time  did  not 
allow  this  additional  processing  to  be  performed  in  the  current 
investigation . 

Finally,  another  aspect  of  the  stalling  process  that 
deserves  some  discussion  is  the  high-frequency  content  that  was 
filtered  out  of  all  the  data  presented  here.  As  mentioned  in 
section  5.3.3,  the  raw  (unfiltered)  Kulite  signals  obtained  at 
70%  speed  and  above  were  dominated  by  a  very  strong  (3-4  psi 
peak-to-peak)  1 /blade  frequency  generated  by  the  passing  shock 
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waves  (see  Fig.  5.38) .  At  60%  speed,  the  entire  span  was 
subsonic  and  the  1/blade  frequency  was  not  as  evident.  It  was 
observed  both  on-line  and  during  data  playback  that  during  the 
transients  to  stall,  the  1 /blade  signals  would  start  modulating 
in  both  frequency  and  amplitude  (perhaps  at  27  Hz) .  The 
modulations  could  be  seen  approximately  1-2  seconds  prior  to 
stall  and  would  grow  as  stall  was  approached.  This  observation 
raises  the  question  of  whether  the  high-frequency  signals  due  to 
shock  passing  (6.9  kHz  at  design  for  the  current  rotors)  might 
provide  a  useful  indicator  of  impending  stall,  particularly  in 
high-speed  machines. 
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6.0  CONCLUSIOHS  AMD  BMCOMMIMDATZONS 


6.1  Sianmry 

Two  single-stage,  transonic  compressor  designs  were  tested 
to  characterize  the  process  leading  up  to  aerodynamic  stall.  The 
two  rotors  had  identical  design  points  (speed,  pressure  ratio, 
mass  flow)  and  were  tested  in  the  same  rig  xinder  the  same 
operating  conditions.  This  allowed  for  comparison  of  the  stall 
inception  behavior  of  the  two  rotors  in  terms  of  the  one 
parameter  that  defined  their  design  difference  -  leading  edge 
sweep  (shock  wave  strength) .  To  quantify  the  stalling  behavior, 
measurements  were  obtained  from  eight  high-response  pressure 
trainsducers  equally  spaced  in  a  plane  0.3  rotor  radii  (2.5 
inches)  upstream  of  the  blade  leading  edge  tips.  Data  were 
obtained  under  both  near-stall  steady  and  transient  compressor 
operating  conditions.  In  addition  to  studying  the  raw  pressure 
traces,  the  signals  were  spatially  emd  teoporally  resolved  using 
Fourier  techniques.  Only  low-frequency  content  was  examined  in 
detail,  and  additional  digital  filtering  was  used  to  enhance 
features  in  the  frequency  range  of  interest  (50%  of  rotor 
frequency) . 

6 . 2  Conclusions 

The  following  conclusions  are  aimed  at  answering  the 
questions  posed  in  section  1.4  of  this  report.  As  is  the  case 
with  most  research  activities,  the  results  of  the  current  work 
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introduced  several  additional  key  issues.  These  are  addressed 
primarily  in  the  recommendations  section. 

(1)  The  transient  stalling  process  of  both  rotors  was 
associated  with  the  growth  of  small  eunplitude  rotating  waves 
(modal  waves)  for  all  speeds  tested.  Initial  wave  an^litude  was 
approximately  0.2  psi  and  did  not  grow  significantly  until  6-10 
rotor  revolutions  prior  to  stall.  The  waves  were  evident 
regardless  of  the  type  of  stall  -  progressive  or  abrupt  -  and  in 
the  presence  of  small  amplitude,  surge-like  oscillations. 

(2)  Compressor  operating  condition  did  not  seem  to  affect 
the  basic  mechanism  of  stall  inception  -  small  amplitude  wave 
growth.  However,  the  length  of  time  the  waves  could  be  detected 
prior  to  compressor  stall  (stall  warning  time)  did  appear  to  be 
related  to  the  shape  of  the  pressure  rise  -  mass  flow 
characteristic  (speedline) . 

(3)  A  positively  sloped  region  of  speedline  characteristic 
appeared  to  be  key  for  detection  of  the  modal  waves.  The  small 
amplitude  rotating  waves  were  not  apparent  in  any  of  the  near¬ 
stall  steady  data  taken  at  speeds  that  exhibited  an  abrupt  stall. 
When  the  compressor  was  operated  transiently  at  the  same  speeds, 
the  rotating  waves  were  detected.  Transient  stall  occurred  at  a 
lower  flow  rate  than  at  the  steady-state  stall  (Figs.  5.35  and 
5.36),  emd  resulted  in  more  of  the  positively  sloped 
characteristic  being  reached.  Additionally,  speedlines  with 
relatively  large  positively  sloped  regions  resulted  in  more 
advanced  stall  warning  times  (Figs.  5.31  and  5.32) . 
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(4)  The  presence  of  shock  waves  in  the  con^ressor  flow 
field  appeared  to  interfere  with  detection  of  the  rotating  waves , 
probably  by  attenuating  the  perturbations.  In  all  cases  where 
the  speedline  shapes  of  the  two  rotors  were  virtually  identical ^ 
the  swept  rotor  (Rotor  6)  exhibited  more  stall  warning  time  them 
the  unswept  design.  The  intent  of  leading  edge  sweep  on  Rotor  6 
was  to  minimize  shock  losses  by  reducing  shock  strength;  the 
reduced  shock  strength*  appeared  to  favor  advanced  detection  of 
the  modal  waves . 

(5)  The  effects  of  speedline  shape  and  shock  waves 
discussed  2d)ove  and  in  sections  5 . 5  amd  5 . 6  are  importamt  to 
sensor  location,  particularly  in  a  multistage  environment.  In  a 
multistage  compressor,  one  would  have  to  be  concerned  with 
individual  stage  characteristics  as  well  as  operating  condition 
(a  multistage  compressor  is  more  likely  to  stall  in  the  front 
stages  during  low-speed  operation  (starting) ;  this  is  not 
necessarily  true  at  high-speed  operation) .  Results  of  the 
current  investigation,  as  well  as  those  from  previous 
investigations  (Gamier  [28]  and  Bonnaure  [20];  see  section  2.4) 
indicate  that  maximum  stall  warning  times  would  require  multiple 
sensor  locations  and/or  advanced  knowledge  of  the  stage  which 
initiates  the  stall. 

(6)  A  considerable  amount  of  digital  data  processing  was 


'Note  that  shock  wave  strength  was  not  characterized  in  the 
current  test  prograuns.  However,  the  swept  rotor  consistently 
exhibited  less  severe  stalls  than  the  unswept  design,  offering 
further  evidence  of  a  weakened  shock  structure  (section  5.6). 
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required  to  obtain  the  stall  warning  times  listed  in  Table  5.1. 
The  processing  included  band-pass  filtering,  spatial  Fourier 
transfozrmation,  data  smoothing,  and  slope  calculation.  Ideally, 
for  use  of  the  modal  waves  as  a  pre-stall  indicator  for  an  active 
engine  control  scheme,  one  would  desire  minimal  processing.  Each 
required  processing  step  would  increase  the  amount  of  stall 
warning  time  needed  to  evaluate  the  signals  and  take  action  to 
avoid  or  suppress  the  impending  stall. 

€ . 3  Recommendations 

Results  presented  in  this  report  have  contributed  to  the 
physical  understanding  of  the  stall  development  process  in  two 
compressor  designs  that  are  typical  of  those  used  in  aircraft  gas 
turbine  engine  applications.  These  results  and  the  potential 
benefits  of  identifying  a  reliable  stall  precursor  (discussed  in 
section  1.1)  warrant  continued  exploration  of  the  stall  inception 
process  in  compressors .  Questions  regarding  the  general 
application  of  the  current  results,  as  well  as  the  robustness  of 
the  small  amplitude  waves,  must  be  addressed  before  any  practical 
applications  can  be  realized.  The  following  recommendations 
outline  a  few  of  these  issues,  with  emphasis  on  the  experimental 
side.  It  should  be  noted,  however,  that  quantitative 
experimental  results  are  essential  if  improvements  in  the 
modelling  of  high-speed  compressor  stability  are  to  be  realized. 

(1)  Stall  inception  investigations  should  continue  on 
various  single-stage  and  multistage  compressor  designs,  with 


66 


emphasis  on  high-speed  machines.  The  stability  theories  appear 
to  predict  quite  well  the  behavior  of  low-speed  machines,  but  are 
in  need  of  quantitative  compressibility  effects  (blade  passage 
choking,  shock  waves,  stage  mismatching,  etc.)  to  accurately 
model  high-speed  compressor  stall  inception. 

(2)  Whenever  possible,  several  axial  locations  for  high- 
response  measurements  should  be  used.  Ideally,  a  measurement  set 
would  be  made  at  the  inlet,  exit,  and  each  compressor  stage  to 
help  quantify  the  effects  mentioned  above.  From  a  practical 
standpoint,  this  would  result  in  very  large  data  storage  problems 
(see  section  4.2)  and/or  very  expensive  analog  recording 
equipment . 

(3)  If  the  equipment  is  available,  storage  of  the 
unfiltered  stall  inception  data  on  magnetic  tape  should  be 
considered  as  the  most  attractive  alternative.  As  was  evidenced 
with  the  current  data  set,  the  advantages  of  magnetic  tape 
storage  -  re-exaunining  raw  signals,  investigating  different 
filter/saunpling  combinations,  etc.  -  proved  to  be  quite  useful, 
and  more  importantly,  provided  additional  information. 

(4)  Exaunination  of  the  high  frequency  content  in  stall 
inception  data  should  be  performed.  A  good  starting  point  would 
be  the  data  set  collected  from  the  current  study.  As  discussed 
in  section  5.7,  the  strong  signal  at  the  blade  passing  frequency 
exhibited  modulations  in  both  frequency  and  aunplitude  as  stall 
was  approached.  The  nature  of  these  modulations  should  be 
studied  to  investigate  their  potential  as  a  stall  warning 
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indication  in  high-speed  compressors.  An  attractive  feature  of 
these  1 /blade  signals  is  that  axial  sensor  location  would  not  be 
as  critical  as  it  appears  to  be  for  detection  of  the  small 
amplitude  modal  waves.  The  strong  1/blade  signal  is  closely  tied 
to  the  passing  shock  waves,  which  are  typically  confined  to  the 
first  (2uid  maybe  second)  compressor  stage. 

(5)  Additional  investigations  should  include  sets  of 
closely  spaced  transducers  to  better  examine  short  length  scale 
phenomena  on  the  stalling  process  (i.e.,  2  or  3  transducers  over 
1  or  2  blade  passages) .  Ideally,  these  measurements  would  be 
obtained  simultaneously  with  the  equally  spaced  measurements,  but 
could  be  performed  separately  provided  adveuiced  provisions  had 
been  made  (i.e.,  holes  drilled  in  the  rig). 

Finally,  it  is  appropriate  to  make  some  comments  that  apply 
specif icly  to  the  Compressor  Aerodynamic  Tvesearch  Laboratory 
(CARL) ,  the  facility  used  in  the  current  study .  The  nature  of 
the  work  conducted  in  the  CARL  regarding  parametric  blade  studies 
makes  it  ideally  suited  for  continued  stall  inception 
investigations .  Transient  data  recording  capabilities  would  have 
to  be  increased  in  order  to  address  the  issues  presented  above, 
2uid  provide  a  complete  picture  of  the  stall  inception  process. 
Tests  conducted  with  different  throttle  closure  rates  and 
distorted  inlet  flow  fields  would  also  be  important. 
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Table  3.1  Facility  Operating  Parameters 


Speed  Range 
Flow  Rzmge 
Rotor  Tip  Distmeter 
Inlet  Total  Pressure 
Steady-State  Pressures 
High  Frequency  Pressures 
Dyneimic  Strain 

Rotor  Tip  Clearamce 


6000  -  21,500  rpm 
20  -  60  Ibm/sec 
14  -  19  inches 
6-15  psia 
160  +  channels 
12  chzmnels 
10  channels  rotating 
20  chamnels  stationary 
8  channels 
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Table  3.2  Rotor/Stator  Design  Features^ 


Design  Speed 
#  of  Stages 
I6Vs 

Tip  Radius^ 

Hub  Radius^ 

Mean  Radius^ 
Rotor/ Stator  Gap 


20,222  rpm 
1 

None 

8.50 

2.65 

5.58 

0.88  tip;  0.29  hxib 


•  llAdoa 

MLA-Cherd 

Caabar* 

Twlafe* 

•allditr* 

Aapaat 

Ratio* 

Rotor  4 

20 

3.95 

12.98 

44.62 

46.86 

2.20 

1.32 

Rotor  6 

20 

4.34 

11.97 

30.10 

41.85 

2.28 

1.25 

Stator 

31 

2.25 

45.30 

14.75 

4.86 

1.78 

1.25 

Notes : 

1.  all  dimensions  in  inches  and  degrees 

2.  at  leading  edge 

3.  at  mean  radius 

4.  Twist  defined  as  tip  stagger  ”  hub  stagger 

5 .  Aspect  Ratio  defined  as  mean  span  divided  by  the  average 
of  chord  at  h\ib,  mean,  and  tip 
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3 . 3  Steady-State  Uncertainty  Analysis 


a.  Mass  Flow  Calculation  Uncertainty 


%Nc 

Parauneter 

Magnitude 

(Ibm/s) 

Uncertainty:  95%  Cl 

Uncertainty:  99%  Cl 

±  Actual 

±  %  Read 

±  Actual 

±  %  Read 

40 

22.611 

0.339 

1.500 

0.474 

2.095 

80 

49.309 

0.264 

0.536 

0.373 

0.756 

100 

61.913 

0.260 

0.420 

0.366 

0.591 

b.  Pressure  Ratio  Calculation  Uncertainty 


%N« 

Parameter 

Uncertainty;  95%  Cl 

Unceirtainty :  99%  Cl 

Magnitude 

±  Actual 

±  %  Read 

±  Actual 

±  %  Read 

40 

1.122 

0.005 

0.404 

0.005 

0.490 

80 

1.580 

0.001 

0.063 

0.002 

0.127 

100 

1.921 

0.001 

0.077 

0.003 

0.161 
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Table  3.4  Kulite  Performance  Characteristics 


Rated  Pressure 
Maximum  Pressure 
Operational  Mode 
Excitation  Voltage 
Maximvim  Excitation 
Sensitivty 

Zero  Pressure  Output 
Temperature  Range 


100  psig 

200  psig 

Differential 

15  VDC 

20  VDC 

1 . 50  mV/psi 

<  ±  5%  Full  Scale 

-40  ®F  to  500  ®F 


Table  3.5  Kulite  Calibration  Data 


Input 

Press . 

(psig) 

Output  Voltage  (V) 

Kulite  # 

1 

2 

3 

4 

5 

6 

7 

8 

m 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2.000 

m 

0.796 

0.783 

0.770 

0.796 

0.851 

0.806 

0.800 

0.791 

1.0 

0.396 

0.376 

0.359 

0.396 

0.458 

0.409 

0.403 

0.386 

m 

0.195 

0.174 

0.155 

0.200 

0.270 

0.208 

0.203 

0.186 

m 

0.075 

0.053 

0.032 

0.080 

0.150 

0.089 

0.085 

0.067 

0.1 

0.032 

0.013 

0.013 

0.040 

0.111 

0.049 

0.044 

0.026 
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Table  4 . 1  Frequency  Rzmge  of  Bemd-paee  Filters 


%  SPEED 

FREQUENCY  BAND-PASS  (HZ) 

60 

60-180 

70 

60-220 

80 

125-230 

90 

130-230 

100 

130-230 

TzJsle  5 . 1  Stall  Warning  Time 

Criterion:  Last  stretch  of  slope  up  to  stall  of  SFC  #1 
between  30-70%  rotor  frequency 


%  SPEED 

ROTOR  4  (UNSWEPT) 

ROTOR  6  (BACKNARD  SWEPT) 

60 

slow,  progressive  stall 

slow,  progressive  stall 

70 

<0.10  sec^ 

slow,  progressive  stall 

80 

<0.10  sec^ 

0 . 2  sec* 

90 

0 . 35  sec* 

>2 . 0  sec* 

100 

0 . 1  sec* 

0.2  sec*'* 

Notes : 

1.  significeuit  period  where  criterion  met  (30-60  rotor 
revs)  well  before  the  stalls  but  separated  by  a 
significant  period  of  random  slope. 

2.  includes  brief  periods  with  slope  outside  criterion. 
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Compressor  Airflow 

e  1.1  Schematic  of  compressor  map 
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Figure  1.2  Compressor  Instabilities  -  surge  and  rotating  stall 
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Figure  1.4  Active  compressor  stabilization 
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Figure  2.4  Spatial  Fourier  coefficient  analysis  of 
modal  waves  [28] 
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Velocity  at  Various  Circumferential  Positions.  Velocity  at  Various  Circumferential  Positions. 


Figure  2.5  Coupling  between  modal  waves  and 
finite  cells  [7] 
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Figure  3.1  Compressor  Aerodynamic  Research  Laboratory  (CARL) 


Figure  3.4  Rotor  4  installed 
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Figure  3.5  Rotor  6  installed 
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Figure  3.7  Close-up  of  Rotor  6  showing  Kulite  location 
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Figure  3.9  Compressor  Instrumentation 
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Figure  3.10  Vane  leading  edge  and  discharge  plane  rake 
instrumentation  [29] 
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Figure  3.12  Kullte  location  and  signal  path 
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%  DESIGN  PRESSURE  RATIO  %  DESIGN  ADIABATIC  EFFICIENCY 
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Figure  4.2 
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Figure  4.3a  Transient  beha>^or  of  an  isolated  airfoil 
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Figure  4.3b  Transient  behavior  of  a  compressor  stage 


Figure  A. 4  Characteristics  of  Butterworth  analog  filter 
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Figure  4.8  Pressure  traces  showing  rotating  waves 
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Figure  4.9  Illustration  of  SFC  #1  analysis  -  90Z  speed  transient 
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frequency  (1/rev) 


Figure  4.10  PSD  of  an  SFC  -  90Z  speed  transient 


eD3OH-f»0>O'HCfT>i(»-o  ffiictomdin'o  — h-cotj 


Time  (rotor  revs) 


b.  digitally  high-pass  filtered 


105 


red  1  !  '  ’  27  Hz 


106 


Figure  5.2  Surge-llke  instability  at  80Z  speed*  Rotor 
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Figure  5.3  Hypothetical  behavior  of  compressor  during 
transient  to  stall  with  planar  oscillations 
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Figure  5.4  Pressure  traces  (no  additional  filtering) 
60Z  speed  transient*  Rotor  4 


l/rev-204  Hz 
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Figure  5.5  Amplitude  and  phase  of  SFC  #1  -  60Z  speed.  Rotor 
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positions  -  60Z  speed.  Rotor 
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Figure  5.7  Pressure  traces  during  70Z  speed  transient.  Rotor  4 
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Figure  5.8  Stall  Inception  at  70Z  speed.  Rotor 
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Figure  5.9  First  and  second  spatial  Fourier  modes  -  70Z  speed.  Rotor 
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Figure  5.12  SFC  #1  of  70Z  speed  transient  data.  Rotor 
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Figure  5«I7  Pressure  traces  during  90X  speed  transient,  Rotor  4 


Time  (rotor  revs) 

Figure  5.18  Stall  inception  at  90Z  speed.  Rotor  4 
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Figure  5.19  SFC  #1  amplitude  and  phase  -  90Z  speed,  Rotor  4 
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Figure  5.20  SFC  #1  amplitude  and  phase  -  90Z  speed.  Rotor  6 
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90Z  speed.  Rotor 
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Figure  5.23  Stall  inception  at  lOOZ  speed.  Rotor  6 
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Figure  5.25  PSD  of  SFC  #1  using  different  lengths  of  transient  data  set  -  lOOZ  speed.  Rotor 
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Figure  5.26  Example  of  smoothed  SFC  #1  phase  data  -  90X  speed 

126 


127 


Figure  5.27  Slope  of  SFC  #1  phase  -  70Z  speed.  Rotor 
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Figure  5.28  Slope  of  SFC  # 
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Figure  5.30  Slope  of  SFC  #1  phase  -  80%  speed.  Rotor 
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Figure  5.31  Slope  of  SFC  #1  phase  -  90%  speed.  Rotor 
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Figure  5.34  Slope  of  SFC  phase  -  100%  speed.  Rotor 
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Figure  5.35  Speedline  characteristics  of  Rotor  4  extended  to 
show  transient  stall  line  observed  during  testing 


Figure  5.36  Speedllne  characteristics  of  Rotor  6  extended  to 
show  transient  stall  line  observed  during  testing 
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Figure  5.37  Effect  of  local  speedline  slope  on  stall  warning  time 
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Figure  5.39  Steady-state  speedline  characteristics  of  both  rotors 
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Appendix 
Fourier  Anelyeie 


Fourier  analysis  plays  an  important  role  in  the  area  of 
digital  (or  discrete)  signal  processing.  Consequently^  in 
addition  to  being  covered  in  many  mathematics  books ^  the  subject 
of  Fourier  analysis  can  be  found  in  virtually  all  signal 
processing  textbooks.  In  this  appendix,  some  basic  Fourier 
concepts  are  reviewed  and  directly  applied  to  the  current  work. 
Detailed  discussions  of  Fourier  techniques  as  they  relate  to 
signal  processing  can  be  found  in  the  textbooks  by  Oppenheim  and 
Schafer  [32]  and  Stearns  and  Hush  [33] . 

Fourier  analysis  allows  for  the  study  of  linear  systems 
based  on  the  representation  of  discrete  and  continuous  time 
functions  in  terms  of  sinusoids  and  harmonics  (multiples) . 
Consider  a  periodic  function,  f (t) ,  composed  of  components  (each 
with  its  own  eunplitude  and  phase)  which  are  sinusoidal  and  occur 
only  at  harmonics  of  the  fundamental  frequency,  (D  Then,  the 
function  can  be  represented  by  a  Fourier  series  in  any  one  of  the 
following  three  forms: 


f(t)  = 


( a  ^coskia  ^  t  +b^sink(a  ^  t ) 

(iCWot+Oj^) 

Ir-l 

K 


-E 


(A.l) 

(A. 2) 


(A. 3) 
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where:  k  »  harmonic  number 

a's,  b's,  A's,  a  'a,  c's  -  Fourier  coefficients^ 
j  ■  imaginary  unit; 
cot  *  2jc/T;  T>period 

The  three  expressions  are  equivalent  as  evidenced  by  the 
relationships  between  the  coefficients  given  in  Table  A.l. 


Table  A.l 

Relationships 

between 

Fourier  Coefficients 

Coefficient 
(k  a  0) 

In  Tenns  of: 

a>.bk 

Ak.  dk 

Ck.C.t 

“t 

b, 

a* 

bk 

Ak  COS  a, 

-At  sin  a. 

c»  +  c-k  (except  oq  =  Ict,) 

j{Ci:  -  C-t) 

Ok 

(ai  +  fci)''" 

Ak 

d» 

2(CkC.k)''^ 

LRelc.)] 

Ck 

■j(a»  -  jbk) 

(t)'~ 

Ck 

<-k 

-j(a»  +  /*») 

c-k 

■  do  ^  0;  Ini.  imaginaty  pan;  Re,  real  pan.  [  33  ] 


In  the  present  work,  the  complex  representation  is  used,  mainly 
for  its  sinq^lif ication  of  the  mathematics .  The  conplex  Fourier 
coefficients  contain  all  of  the  spectra  information  (distribution 
of  amplitude,  phase,  power,  or  energy  over  frequency)  of  the 
signal . 

Now,  consider  the  case  of  the  present  investigation 
indicated  in  Fig.  3.12.  The  eight  transducers  are  equally  spaced 
around  the  conqpressor  annulus  (45  degrees  apart) ,  which  implies 

^Also,  least~squares  coefficients,  as  Fourier  series  are  an 
application  of  least-squares  approximations. 
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the  functions  comprising  the  Fourier  representation  are 
orthogonal  (i.e.,  the  exponentials  in  Eqn.  A.3) .  In  order  to 
determine  the  presence  of  modal  waves /  the  Spatial  Fourier 
Coefficients  (SFCs)  must  be  computed  at  each  discrete  time.  The 
spectra  of  the  seunpled  signals  can  then  be  tracked  over  the  total 
sample  time.  The  SFCs  can  be  formulated  using  the  solution  for 
least-squares  coefficients  (Stearns  and  Hush  [33])  to  obtain: 

N-l  -d2itnk 

n-O 

where:  p(0„,t)  unsteady  pressure  measurement  at  a  particular 

circumferential  location  and  time 
N  »  number  of  spatial  measurements;  N«8  here 
k  =  harmonic  *  ±1,  ±2,  ±3* 

Equation  (A.  4)  represents  the  discrete  Fourier  transform 
(DFT)  of  a  periodic  signal.  Since  the  measurements  are  equally 
spaced  and  the  measured  pressures  are  real^  the  special  symmetry 
properties  of  the  DFT  can  be  used;  namely,  c,,  and  c_,j  are  complex 
conjugates.  Thus,  only  c_i,  c.j,  and  c.j  need  be  computed  (or  the 
positive  harmonics)  to  determine  the  desired  amplitude  and  phase 
information.  It  is  easily  seen  from  Table  A.l  that  the  c_,, 
contain  the  phase  and  one-half  the  aunplitude  information  of  each 
harmonic . 

In  addition  to  amplitude  and  phase  ,  the  SFCs  can  be  used  to 


^Fourier  analysis  states  that  in  order  to  confute  K 
harmonics,  the  number  of  measurements  N  must  be  2R+1 .  Thus,  with 
8  measurements,  only  the  first  3  harmonics  can  be  resolved. 
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compute  the  power  spectral  densities  (PSDs)  of  each  Fourier  mode. 
The  PSD  of  each  SFC  is  proportional  to  the  square  of  the 
amplitude  of  the  SFC,  |Ck|^.  As  discussed  by  Paduano  [34],  for 
applications  involving  analysis  of  modal  waves,  care  must  be 
talcen  in  interpreting  the  PSD  plots.  Because  the  SFCs  are 
complex  and  represent  a  spatial  wave,  a  frequency  on  a  PSD  plot 
can  represent  a  rotation  frequency  of  the  wave  (which  is  of 
interest  here)  or  a  frequency  of  oscillation  in  wave  amplitude. 

A  PSD  which  is  symmetric  aOsout  zero  frequency  will  represent  the 
latter  behavior.  Thus,  only  PSD  peaks  which  do  not  exist  at 
identical  positive  and  negative  frequencies  indicate  the  modal 
waves.  As  suggested  by  Paduauio,  overlapping  the  positive  and 
negative  frequency  plots  provides  a  convenient  representation  of 
the  PSD. 
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